T/Hb-flififi  195  IONIC  MECHANISMS  OF  SOOT  FORMATION  IN  FLAMES(U) 
flEROCHEM  RESEARCH  LABS  INC  PRINCETON  NJ 


UNCLASSIFIED 


H  F  CALCOTE  ET  AL 
F49620-85-C-0150 


JUN  87  AFOSR-TR-87-1197 


(me  FlfF  mpy 
AD-A186  195  _ 


|J  1 1  i  J  |  ^ 


«foat  s*cuaitv  classy  cation 
.class  tf  led 


UMCNTATION  PAGE 


1b.  RESTRICTIVE  MARKINGS 

None 


:y«iTY  Classification  a 


»  Ot Classification  oownsff 


4  F|F*OKMINQ  ORGANIZATION 


il 


3  O'STRl  BuT  ION/ A  V  Ai  LAS  I  LI  TV  OF  REPORT 

Distribution  unlimited;  approved  for 
public  release 


(X>  0 


S.  MONITOR! 


*A  NAM*  OF  FCRFORMlNG  ORGANIZATION 

AeroChem  Research  Labs.,  Inc. 


•c.  AOORISS  C|I>  jiui  «<  /If  CMi 

P.0.  Box  12 
Princeton,  NJ  08542 


Sf.  NAM*  OF  F  gNOiNG  SPONSORING 
ORGANIZATION 


7«.  NAM*  OF  MONlTORlNO  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 


Tb_AOOR*SS  iCifyi 

Bering  AF1 


IS  (Cit*. IFF  «m4  ZIP  CoOti 

ng  AFB  DC  20332-6448 


OFFlCt  SYMBOL  Is  FROCURCMINT  INSTRUMSNT  iQCNTiFiCATION  NUMSCR 


Air  Force  Office  of  Sci.  Res.  I  AFOSR/NA 


F49620-83-C-0150 


11  TiTlI  Inciue*  Security  C^m*hcoitoni 

Ionic  Mechanisms  of  Soot  Formation  in  Flames 


U  FIRSONAL  AijTHORlSl 

H.F.  Calcote  and  D.G.  Ke i 1 


10.  SOURCE  OF  FUNDING  NOS 

PROGRAM 

FROJSCT 

TASK 

WORK  UNIT 

«l*ment  no. 

NO. 

NO. 

NO 

61102F 

2308 

A2 

13a  t*m  of  report 

13A  Tim*  COVCRSO 

14.  OAT*  OF  REFORT  (Yr.  Mo..  Dot! 

|  Final 

mom 

June  1987 

1  1 1  7 

COSATi  < 

!  LjTi 

:  LQ 

1  GAOUF 

IS  SUBJECT  TIRMS  i  C  on  tin  wc  of  iyikim  if  noeomary  and  tOanttfy  by  Woe*  iwntri 

Soot  Formation;  Ionic  Mechanisms;  Flame  Measurements;  Mass 
Spectrometry;  Langmuir  Probes;  Thermocouple  Probes; 


19.  ABSTRACT  /Continue  on  if  neetmory  end  identify  by  btoeb  number  f 

Experimental  measurements  have  been  made  and  interpreted  in  acetylene/oxygen  and 
benzene/oxygen/argon  flames  at  2.7  kPa,  and  an  unburned  flow  velocity  of  50  cm/s 
with  the  objective  of  evaluating  the  ionic  mechanism  of  soot  nucleation.  This 
mechanism  postulates  that  chemi-ions  are  the  precursors  of  soot  and  that  the 
initial  reactions  in  the  soot  nucleation  process  are  ion-molecule  reactions  in 
which  molecular  ions  continually  increase  in  size  until  they  are  neutralized  by 
ion  recombination.  Total  ion  profiles  were  determined  by  Langmuir  probe;  indiv¬ 
idual  ion  profiles  were  determined  by  molecular  ion  sampling  mass  spectrometry  up 
to  about  mass  600;  temperature  profiles  were  determined  by  radiation  corrected 
thermocouples  It  is  demonstrated  that  the  ion  concentration  peaks  ahead  of  the 
appearance  of  soot;  the  ion  concentration  is  greater  than  the  concentration  of 
soot  particles;  and  ions  decay  as  soot  is  produced.  In  the  acetylene/oxygen 
<t>  *  3  flame,  the  ion-molecule  reaction  rates  are  measured  and  compared  with  other 


(over) 


30  O'STRl  SUTlON/A  V  Ai LAOI  LI T Y  OF  ABSTRACT 
JNCLASSlFiEO/UNLiMITEO  X  SAM*  AS  RFT.  C  OTIC  US* RS  □ 


33a  NAM*  OF  R*SFONSi*L(  iNOiVIOUAL 

Julian  M.  Tishkoff 


OO  FORM  1473,  83  APR  •  OlTION  OF  1  JAR 


31.  ABSTRACT  SSCURITV  CLASSIFICATION 

Unclassified 


I 


33b.  T* LCFHONC  NUMSIR 
IfncluOt  Art*  Co4t) 

(202)767-4935 


•  OlTION  OF  1  JAN  T3  IS  OSSOL*T*.A  _  - 

87  9  24  206 


133c.  OFFICE  SYMBOL 

AFOSR/NA  _ 


Unclassified _ 

SSCURITV  CLASSIFICATION  OF  THIS  FACE 


UNCLASSIFIED 


»,*1 


,»»'i 

'»!] 


u 

A 


**  j 

'*!J 


SCCUMiTv  classification  of  this  faqc 


BLOCK  18  (Continued) 


Ion-Molecule  Reaction  Rates 


BLOCK  19  (Continued) 

measured  and  calculated  rates.  It  is  demonstrated  that  these  rates  are  rapid  at  flame 
temperatures.  Some  major  differences  were  found  in  the  features  of  benzene  and 
acetylene  ion  profiles  that  remain  to  be  explained.  The  experiments  are  strongly 
supportive  of  the  ionic  mechanism  of  soot  formation  in  flames. 


UNCLASSIFIED 


SCCUMITV  CLASSIFICATION  OF  THIS  FAOS 


&FOSR.TR.  8  7-1197 

TP-465 

TABLE  OF  CONTENTS 


DD1473  » 

I .  INTRODUCTION  1 

II.  APPARATUS  AND  EXPERIMENTAL  TECHNIQUES  2 

A.  Apparatus  2 

B.  MIT  Burner  2 

C.  Temperature  Measurements  2 

D.  Mass  Spectrometer  Mass  Calibration  8 

E.  Calibration  of  Mass  Spectrometer  Sensitivity  8 

III.  RESULTS  AND  DISCUSSION  9 

A.  Temperature  Effects  9 

B.  Thermodynamics  of  Soot  Formation  10 

C.  Reaction  Heating  of  Soot  Particles  11 

D.  Acetylene-Oxygen  Flames  IS 

1.  Total  Ion  Concentrations  16 

2.  Ion  Spectra  and  Mechanistic  Implications  16 

3.  Ion  and  Neutral  Species  Profiles  and  Reaction 

Rate  Implications  18 

E.  Benzene-Oxygen  Flames  22 

1.  Total  Ion  Concentrations  23 

2.  Calibration  of  Mass  Spectrometer  Sensitivity  25 

3.  Mass  Spectra  25 

IV.  PUBLICATIONS  25 

V.  PERSONNEL  26 

VI.  TECHNICAL  INTERACTIONS  27 

VII.  INVENTIONS  AND  PATENT  DISCLOSURES  29 

VIII.  REFERENCES  29 


1  Suggested  Structures  of  Some  of  the  Observed  Ions 


34 


II  Occurrence  of  Carbon  and  Hydrogen  Atoms  in  Flaae  Ions  36 

III  Rate  Constants  for  Reactions  of  Neutrals  with  Ions  to  Produce 

Larger  Ions  37 

LIST  OF  FIGURES 

F i qure  Paoe 

1  THERMOCOUPLE  HEATING  CURVES  IN  VACUUM  AND  IN  0  ■  2.25  FLAMES  38 

2  EXPERIMENTAL  TEMPERATURE  PROFILES  FOR  ACETYLENE-OXYGEN 

FLAMES  FOR  DIFFERENT  EQUIVALENCE  RATIOS  38 

3  COMPARISON  OF  TEMPERATURE  PROFILES  IN  THE  MIT  DESIGN  COPPER 

BURNER  AND  THE  AEROCHEM  STAINLESS  STEEL  MULTITUBE  BURNER  39 

4  TEMPERATURE  PROFILES  AS  FUNCTION  OF  DISTANCE  ABOVE  BURNER 

NORMALIZED  BY  BURNER  DIAMETER  39 

5  MASS  CALIBRATION  FUNCTION  40 

6  COMPARISON  OF  TOTAL  ION  CONCENTRATION  PROFILE  FROM 
LANGMUIR  PROBE  MEASUREMENTS  WITH  PRESENT  MASS 

SPECTROMETER  MEASUREMENTS  41 

7  COMPARISON  OF  ADIABATIC  AND  MEASURED  TEMPERATURES  AS  A 
FUNCTION  OF  SOOT  THRESHOLD t  0« ,  FOR  ATMOSPHERIC 

TOLUENE/O, /N,  FLAMES  42 

8  EQUILIBRIUM  COMPOSITIONS  CALCULATED  AT  THE  MEASURED 

SOOT  THRESHOLD  FOR  A  SERIES  OF  FUELS  43 

9  COMPARISON  OF  EQUILIBRIUM  CALCULATIONS  OF  CHARGED  SOOT 

(POSITIVE  ION)  CONCENTRATIONS  WITH  MEASURED  CONCENTRATIONS  43 

10  INCREASE  IN  ION  MASS  WITH  DISTANCE  ABOVE  BURNER 

DEDUCED  FROM  HOMANN  AND  STROEFER'S  MEASUREMENTS  44 

11  INDIVIDUAL  ION  PROFILES  45 


TP-465 


*,* 

i\r 


>>»< 

w 

w 

fV> 

#*4* 

it  j 


t*  a 

$ 

I 

•?< 

J 

yi 

4 

a 

sa 


a 

*£| 


F i auy 


Plflt 


CONCENTRATION  PROFILES  OF  SELECTED  NEUTRAL  AND 
IONIC  SPECIES 

COMPARISON  OF  MEASURED  (POINTS)  WITH  CALCULATED  (LINES) 
DIFFUSION  COEFFICIENTS 

LANGMUIR  PROBE  CURRENT  ANO  ION  CONCENTRATION  PROFILES 

TEMPERATURE  PROFILES  IN  0  «  1.8  FLAME 

LANGMUIR  PROBE  CURRENT  PROFILE 

HIGH  PASS  FILTER  MASS  SPECTROMETER  PROFILES 

HIGH  PASS  FILTER  MASS  SPECTROMETER  PROFILES 

MASS  SPECTROMETER  SENSITIVITY  CALIBRATION  FUNCTION 

INDIVIDUAL  IONS  IN  0  -  1.8  FLAME 

INDIVIDUAL  IONS  IN  0  -  1.8  FLAME 

INDIVIDUAL  IONS  IN  0  «  1.8  FLAME 


INDIVIDUAL  IONS  IN  0  -  1.8  FLAME 


INDIVIDUAL  IONS  IN  0  ■  1.8  FLAME 


LOCATION  OF  CURRENT  MAXIMA  IN  ION  PROFILES 


MAXIMUM  CURRENTS  IN  ION  PROFILES 


INDIVIDUAL  ION  PROFILES  IN  0  >  2.0  FLAME 


Accedo'i  '-‘.jr 

NTIS 

DT!C  7  •  R  [) 

Ui.2ijr.o  .  ;  !  ;  j 

JiictitiC.  1 


By . 

Dist.-it,  / 


1 

— ] 


a 

I*?, 


y< 


Avj'i  .i”  i or 


k-/ 


7*1 *i  V  •  ‘ 

*,  4.4' 


For  a  number  of  years  tha  personnel  at  AaroChaa  hava  baan  protagonists  of 
tha  ionic  aechanisa  of  soot  foraation  in  flaaas.  In  this  theory,  it  is  as- 
suaad  that  cheai-ions  grow  vary  rapidly  through  a  sarias  of  ion-aolacula  ra ac¬ 
tions  to  produca  larga  ions  which  ara  nautralizad  to  produca  aithar  vary  saall 
incipiant  soot  particlas,  or  vary  larga  nautral  aolaculas  which  rapidly  grow 
to  fora  tha  incipiant  soot  particlas.  Ua  hava,  in  fact,  previously  deaon- 
stratad  that  thara  is  no  sharp  daaarcation  between  larga  aolaculas  and  saall 
soot  particlas. * »* 

After  auch  study  by  aany  people,  an  understanding  of  tha  soot  nucleation 
process  still  represents  one  of  tha  aore  significant  problaas  challenging  tha 
coabustion  coaaunity.  Much  of  tha  reason  for  this  is  that  tha  necessary  data 
required  to  resolve  tha  aany  questions  about  tha  initial  stages  of  soot  foraa¬ 
tion — where  tha  major  questions  apply — ara  vary  difficult  to  obtain  experimen¬ 
tally}  global  aeasureaents  of  soot  particlas  do  not  give  tha  necessary  infor- 
aation.  Tha  objective  of  our  program  has  baan  to  make  the  necessary  detailed 
ion  (both  total  and  individual)  concentration  profile  aeasureaents  in  tha  same 
flaaas  in  which  Howard  and  associates  at  MIT1**  and  others,  e.g.,  Del f au  and 
associates1**  in  Franca,  and  Hoaann  and  associates  in  Germany****1  ara  making 
complementary  aeasureaents.  Our  aeasureaents  included:  mass  spectrometer 
aeasureaents  for  individual  ion  concentrations!  Langmuir  probe  aeasureaents 
for  total  ion  concentrations}  and  temperature  compensated  thermocouple  aea¬ 
sureaents  for  temperature  profiles.  To  assure  that  the  flames  we  study  are 
exactly  comparable  to  those  used  by  Howard  and  associates  at  MIT,  we  have 
duplicated  their  burner. 

The  work  statement  for  this  prograa  contained  two  phases.  Phase  I  exper¬ 
iments  were  to  be  done  in  well  studied  flames  of  benzene-oxygen  and  acetylene- 
oxygen  and  they  included  generation  of  temperature  profiles,  Langmuir  probe 
curves  of  total  ion  concentration,  and  mass  spectrometer  profiles  of  individ¬ 
ual  ions.  Phase  II  involved  model  development  of  the  individual  processes 
involved  in  incipient  soot  foraation,  e.g.,  rates  of  ion  and  particle  foraa¬ 
tion,  rates  of  production  of  large  molecular  ions,  the  temperature  of  the 
growing  particle  and  its  effect  on  thermal  ionization,  thermal  ionization  of 
large  molecules,  and  ion  molecule  rate  coefficients.  The  quantitative  models 
developed  in  this  program  for  each  of  the  steps  will  fora  the  basis  for  a 
detailed  overall  quantitative  model  of  the  total  process  of  soot  foraation. 
Considerably  aore  tiae  was  consumed  in  the  calibration  of  the  mass  scale  of 
the  aass  spectrometer  than  anticipated  but  this  was  necessary  in  order  to 
treat  the  problea  quantitatively.  This  was  done  at  the  expense  of  other  ex¬ 
perimental  tasks. 

In  this  report  we  summarize  the  results  on  this  prograa.  Results  which 
have  been  published  are  only  briefly  summarized}  other  results  are  covered  in 
aore  detail. 
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IX.  APPARATUS  AND  EXPERIMENTAL  TECHNIQUES 


A.  APPARATUS 

The  apparatus  used  in  the  aass  spectrometric  studies  consists  of  a  multi- 
tubular  burner  supported  on  a  vertical  translation  stage  which  feeds  through  a 
vacuum  seal  into  a  cylindrical  flame  chamber  pumped  by  a  large  (140  L  s* 1 ) 
mechanical  pump.  The  top  of  the  chamber  is  fitted  with  a  water-cooled  plate 
supporting  a  metal  sampling  cone  which  permits  the  sampling  of  flame  gases 
through  an  (0.013-0.025  cm)  orifice  into  the  first  of  two  differentially 
pumped  vacuum  chambers.  This  chamber,  maintained  at  less  than  0.01  Pa  pres¬ 
sure,  contains  a  series  of  electrostatic  lenses  which  focus  and  direct  the 
flame  ion  beam  into  the  second  vacuum  chamber  containing  a  quadrupole  mass 
filter,  maintained  at  pressures  below  0.001  Pa  by  a  IS  cm  (6  in.)  diffusion 
pump.  The  quadrupole  used  in  this  work  was  constructed  at  AeroChem  about  20 
years  ago.  In  this  study  it  was  powered  by  a  350  kHz  rf  supply  providing  a 
useful  aass  range  up  to  nearly  600  amu.  Ions  which  pass  through  the  mass 
filter  strike  the  cathode  of  an  electron  multiplier  and  the  amplified  current 
is  recorded  on  an  XY  recorder  as  a  function  of  mass  (fixed  burner  position)  or 
as  a  function  of  burner  position  (fixed  ion  mass).  A  digitizer  pad  was  used 
to  transfer  mass  spectrometer  data  from  recorder  charts  to  a  computer  for 
processing. 


B.  HIT  BURNER 

The  copper  burner  used  by  Bittner  and  Howard11  at  HIT  was  duplicated  from 
designs  furnished  by  Howard14  so  we  could  compare  our  ion  profiles  with  their 
profiles  of  neutral  species  in  the  same  flames.  It  consists  of  three  major 
pieces.  The  burner  surface  is  a  7.1  cm  diameter,  1.27  cm  thick  copper  plate 
with  about  600  1  am  holes  drilled  in  a  hexagonal  close-packed  pattern  such 
that  the  center  to  center  distance  is  0.254  cm.  This  plate  is  press-fit  into 
a  thin-walled  (0.3  cm  thick)  support  cylinder,  which  is  threaded  onto  the 
bottom  mixing  chamber  and  water-cooling  chamber.  Therefore  the  burner  surface 
is  cooled  by  conduction  down  the  copper  support  cylinder  to  the  water  jacket 
at  the  burner  base.  This  burner  operates  with  a  surface  temperature  (*  425  t 
25  K  at  burner  edge)  which  is  greater  than  our  directly  cooled  stainless  steel 
burner.14  We  had  considerable  difficulties  getting  this  burner  to  operate  in 
our  system)  the  flame  showed  a  slow  periodic  oscillation.  We  finally  traced 
this  problem  to  a  gas  mixing  chamber  and  eliminated  the  difficulty  by  simply 
changing  the  manner  in  which  the  two  gases  were  brought  into  the  chamber.  No 
similar  problem  was  observed  with  our  stainless  steel  burners  which  have  larg¬ 
er  pressure  drops  across  the  burner. 


C.  TEHPERATURE  MEASUREMENTS 

Accurate  temperatures  are  required  to  interpret  the  various  data  to  be 
analyzed  and  they  are  generally  not  available  in  the  flames  in  which  we  were 
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interested.  Measurement  of  flame  temperatures  in  sooting  flames  is  a  compli¬ 
cated  problem  and  is  thus  often  avoided  or  done  without  correcting  for  radia¬ 
tion  losses.  Me  thus  report  our  procedures  here  in  some  detail;  we  plan  to 
publish  this  work  in  the  future. 

For  our  primary  measurements  we  chose  to  use  a  thermocouple  technique  in 
which  radiation  losses  are  compensated  for  by  electrical  heating  of  the  ther¬ 
mocouple,  but  because  of  this  technique's  complexity,  we  also  used  it  to 
calibrate  the  radiation  corrections  of  coated  thermocouples  which  were  not 
electrically  heated.  As  a  check,  these  radiation  corrections  were  also  calcu¬ 
lated  by  the  usual  techniques.  Fine  wire  S-type  and  B-type  thermocouples  1  to 
2  cm  long  with  diameters  from  50  to  125  pm  were  fused  onto  250  pm  wires  of  the 
same  materials  as  saddle  supports.  The  support  wires  fed  through  an  alumina 
tube  which  was  in  turn  inserted  into  a  stainless  steel  0.64  cm  o.d.  tube, 
vacuum-sealed  with  epoxy  cement,  and  brought  through  the  low  pressure  burner 
housing  wall  through  a  sliding  0-ring  seal. 

In  nonsooting  flames  catalytic  heating  of  the  thermocouple  was  avoided  by 
coating  the  thermocouple  with  Be0/Y203.*7  The  coating  was  fused  at  about  1870 
K  in  a  Meeker  flame.  In  sooting  flames  the  soot  which  rapidly  formed  on  the 
probe  was  found  to  be  a  good  noncatalytic  coating.  In  some  cases  the  thermo¬ 
couple  was  intentionally  soot -coated  in  a  rich  flame  and  rapidly  used  in  a 
nonsooting  flame  before  the  coating  burned  off.  The  diameters  of  coated  wires 
were  measured  before  and  after  the  experiments  using  a  microscope  to  determine 
the  coating  thickness  which  was  typically  about  10-15  pm. 

In  the  electrically  heated  thermocouple  method,  the  thermocouple  tempera¬ 
ture  is  measured  in  a  vacuum  (S  3  x  10- a  Pa)  as  a  function  of  heating  current. 
This  is  done  by  resistively  heating  the  wire  with  a  3  kHz  a.c.  current  and 
simultaneously  measuring  the  d.c.  thermal  emf  of  the  thermocouple  junction.  In 
a  good  vacuum  the  convection  losses  are  small  or  negligible  so  that  the  elec¬ 
trical  power  loss  in  the  wire  is  equal  to  the  radiation  losses.  A  plot  of 
junction  temperature  vs.  heating  current  is  thus  generated.  The  experiment  is 
then  repeated  in  the  flame  producing  another  (different)  plot  of  thermocouple 
temperature  vs.  heating  current.  If  the  two  curves  (one  obtained  in  vacuum 
and  one  in  the  flame)  are  plotted  on  the  same  coordinate  system  they  will 
cross  at  some  point.  If  the  radiative  properties  of  the  thermocouple  (6  and 
d« ,  where  €  *  the  wire  emissivity  and  d«  =  the  thermocouple  wire  diameter)  do 
not  change  between  the  two  experiments,  the  temperature  at  the  crossing  point 
is  the  true  flame  temperature.  At  this  point  there  are  no  net  losses  from  the 
thermocouple.  To  assure  that  the  radiation  properties  of  the  thermocouple  do 
not  change  during  each  experiment,  the  temperature-current  curve  is  repeated 
in  vacuum  after  each  flame  experiment. 

The  above  procedure  is  tedious,  and  due  to  the  fragile  nature  of  the  wire 
coating  it  was  necessary  to  produce  heating  curves  rapidly  before  the  surface 
properties  changed.  Many  thermocouples  were  lost,  especially  at  the  highest 
temperatures  so  that  extrapolation  procedures  were  developed. 
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Under  vacuum  conditions,  in  the  absence  of  convection  and  conduction 
losses  at  the  thermocouple  junction,  the  energy  balance  involves  only  elec¬ 
trical  heating  and  radiative  losses  to  the  surroundings.  Equating  the  rates 
per  unit  wire  length! 


a€ndt  (T»*  -  T**) 


ltd* 

w 


1* 


where  a  is  the  Stefan-Boltzmann  constant,  d,  is  the  total  wire  plus  coating 
diameter,  d»  is  the  diameter  of  the  wire  itself,  p  is  the  resistivity  of  the 
junction  area  and  1  is  the  heating  current.  The  wire  temperature,  T» ,  is  much 
greater  than  the  temperature  of  the  surroundings,  T»  so  that 


«€ndr  T»  4  -  1* 


The  current  required  to  heat  the  thermocouple  to  a  temperature  T»  is  then: 


For  platinum  and  rhodium  alloys,  both  €  and  p  show  similar  temperature  depen¬ 
dencies  in  the  range  of  interest,  1000-2000  K.  Therefore  the  overall  current 
to  temperature  proportionality  is  I  «  T» * . 

The  calculated  slopes  dl/dT*  should  provide  a  measure  of  (Cd,)0**.  The 
experimental  values  were  grouped  by  the  thermocouple  visual  appearance  and 
dl/dT*  normalized  (to  a  new,  shiny  thermocouple).  Slopes  are  given  here: 

_ Appearance _  <dI/dTM»)/ (dl/dT, 


Bare,  shiny 
Bare,  not  shiny 
Bare,  rough 

Dark  coating,  smooth  -  same  diameter 
Dark  coating,  rough  -  same  diameter 
Clean  Be0/Ys03  -  thick 
Dark  Be0/Y*03  -  thick 


1.00 

1.43 

1.57  ±  0.08 
1.80  ±  0.10 
1.86  ♦  0.05 
1.96 
2.36 


This  provides  a  rationale  for  extrapolation  of  the  vacuum  heating  curves. 

Over  20  heating  curves  in  vacuum  were  made  for  variously  treated  thermocoup¬ 
les.  The  current  required  to  heat  the  thermocouples  to  any  fixed  temperature 
varied  by  nearly  a  factor  of  three,  yet  a  linear  least  squares  fit  of  T«*  vs. 
I,  always  gave  a  correlation  coefficient  r*  >  0.99.  The  T»*  vs.  X  fitting 
routine  was  thus  used  to  extrapolate  to  temperatures  not  directly  accessible. 
The  heating  curves  in  the  flame  environment  do  not  lend  themselves  to  such 
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simple  closed  form  analysis  since  the  convective  heat  transfer  to  the  wire  has 
a  different  temperature  dependence  from  the  radiation. 

Three  sets  of  heating  curves  are  shown  in  Fig.  1  for  two  different  ther¬ 
mocouples  (circles  and  triangles)  in  a  0  *  2.25  flame,  21  mm  above  the  burner 
and  for  one  thermocouple  (squares)  1.5  mm  above  the  burner  in  the  same  flame. 
Thermocouple  A  was  more  heavily  coated  with  carbon  or  soot  for  these  measure¬ 
ments.  Both  thermocouples  were  B-type  using  51  pm  wires.  The  open  symbols 
represent  the  vacuum  calibration  data  while  the  solid  symbols  represent  the 
heating  curves  in  the  flame.  The  curves  through  the  open  symbols  are  (linear) 
least  squares  fits  to  the  data  above  1000  K  to  the  form:  I  =  a  ♦  bT*2.  (The 
line  through  the  open  circles  represents  a  linear  extrapolation  of  the  four 
data  points  above  1300  K.)  Either  extrapolation  method  gives  about  the  sane 
result  at  the  intersection  of  the  vacuum  and  flame  heating  curves.  The  flame 
data  (solid  symbols)  exhibit  curvature  at  low  currents.  However,  as  the  wire 
temperature  increases,  the  T«  heat  loss  term  becomes  more  important  and  the 
curvature  decreases.  Simulated  heating  curves  calculated  here  have  shown  that 
the  approach  to  the  flame  temperature  is  well  approximated  by  a  linear  temper¬ 
ature  vs.  a.c.  current  relationship.  Ue  feel  that  a  linear  extrapolation  of 
the  flame  data  is  more  realistic,  although  it  is  recognized  that  uncertainties 
in  extended  extrapolation  can  be  fairly  large.  Figure  1  supports  the  extrapo¬ 
lation  method.  The  two  thermocouples  used  at  21  mm  have  significantly  differ¬ 
ent  heating  curves,  yet  the  extrapolated  flame  temperatures  are  the  same  for 
both  cases. 

For  thermocouples  which  are  not  electrically  heated  (1=0),  in  the  ab¬ 
sence  of  catalytic  heating  and  wire  conduction  losses  (negligible  with  the 
long  («  1  cm)  fine  wires  used  here),  the  heat  balance  equation  gives  the  tem¬ 
perature  correction  for  radiation  for  cylindrical  thermocouples, 


«€ 

AT  s  T,  x  -  T  *  —  T« 

■  h  * 


where  T,  t  is  the  true  local  flame  temperature,  T»  is  the  thermocouple  junction 
temperature,  and  h  is  the  heat  transfer  coefficient  from  the  flame  gases  to 
the  wire.  The  temperature  of  the  surroundings  to  which  the  thermocouple 
radiates  is  assumed  to  be  small  relative  to  T*..  The  difference  between  the 
flame  temperature  from  the  intersection  of  the  heating  curves  and  the  thermo¬ 
couple  reading  with  no  heating  current  is  the  radiation  correction,  AT.  The 
parameter  h  can  be  expressed  as 


kr  Nu 


where  k,  is  the  thermal  conductivity  of  the  flame  gases  and  Nu  is  the  Nusselt 
number.  Thus, 
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€dt  ffT,  « 


Under  conditions  appropriate  to  low  Reynolds  number,  Re,  as  in  these  measure¬ 
ments,  Nu  has  been  approximated  as1*  0.8(Re>°-**  and  as1*  (0.42  Pr°- *  ♦  0.57 
pro- 33  Re°- 5)  where  the  Prandtl  number  Pr  is  roughly  unity.  In  either  case, 
Nu  is  only  mildly  dependent  on  the  thermocouple  dimensions  (Re  «  dT ) . 

Thus,  for  the  same  position  in  a  flame  we  estimate  that  the  thermocouple 
response,  in  terms  of  AT/T»4,  is  proportional  to  (dT€).  The  vacuum  heating 
curves  were  shown  above  to  scale  as  1  «  (dt€)°-*T»*.  Therefore  the  radiation 
correction  should  be  closely  related  to  the  square,  (dI/dT*>*.  For  the  two 
thermocouples  used  at  21  mm  (Fig.  1)  the  ratio  of  the  correction  term  AT/T. * 
is  about  0.7  while  the  ratio  of  (dl/dT**)*  for  the  vacuum  calibrations  is 
about  0.6,  consistent  with  the  predictions. 

In  addition  to  the  above  procedures  the  heat  transfer  coefficient  was 
calculated  based  on  Kaskan’s  treatment1*  which  assumes  Nu  *  0.8(Re)<>*  **. 
Estimates  of  the  flame  composition  were  obtained  from  equilibrium  calcula¬ 
tions.  Viscosity  was  calculated  with  Wilkes  mixture  rule  and  thermal  conduc¬ 
tivity  was  calculated  with  Uissiljewa's  mixture  rule  as  described  in  Ref.  20. 
For  the  three  calibrations  in  Fig.  1  the  experimental  values  of  AT  and  the 
calculated  values  for  h  were  consistent  with  a  reasonable  emissivity,  €  =  0.8. 
Due  to  uncertainties  in  the  values  of  h,  we  do  not  claim  to  have  measured  € 
but  an  elevated  emissivity  above  that  for  a  shiny  thermocouple  is  consistent 
both  with  visual  observations  of  carbonaceous  coatings  formation  in  these 
nonsooting  flames  and  with  the  vacuum  calibrations. 


Determination  of  flame  temperatures  with  thermocouple  heating  curves  is 
both  time  consuming  and  tedious.  However,  one  can  make  use  of  any  single 
measurement  to  determine  €/h  for  a  thermocouple  in  a  particular  flame  as  de¬ 
scribed  above.  Since 


€ 

h 


kr  Nu 


one  can  estimate  the  variation  in  €/h  and  thence  the  radiation  correction,  AT. 
The  Nusselt  number  is  assumed  to  only  weakly  depend  on  the  local  flame  proper' 
ties.  For  a  given  thermocouple  with  a  fixed  surface,  €dT  is  constant.  Note 
that  dI/dT»*  is  a  measure  of  this  for  any  specific  thermocouple  and  can  be 
used  to  check  the  thermocouple  properties.  The  thermal  conductivity  kr  is 
composition  and  temperature  dependent.  However,  the  major  product  gases  in 
0  *  2.5  to  3.5  flames  have  similar  thermal  conductivities,  so  the  calculated 
temperature  dependence  of  kr  of  about  T1  is  the  major  correction  factor. 
Electrically  heated  thermocouple  temperature  measurements  were  made  near  the 
location  of  the  maximum  temperature,  at  T»*» ,  in  the  0  *  3.0  flame  (at  1.0  cm 
from  burner).  The  values  of  1883  and  1923  K  were  determined  with  two  differ¬ 
ent  thermocouples.  The  differences  between  the  uncorrected  thermocouple  read 
ing  in  the  flame  and  the  flame  temperature  was  on  the  order  of  270-300  K.  Tw< 
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values  of  €/h  were  then  calculated  and  the  mean  value  was  used  to  correct  all 
the  flame  profiles  for  2.5  £  0  £  3.5,  correcting  for  the  temperature  depen¬ 
dence  of  k,  scaled  as  T1.  With  correction  factors  on  the  order  of  200  to  400 
K,  the  inclusion  of  temperature  dependent  k,  normally  changed  the  correction 
by  less  than  40  K.  The  effect  of  flame  composition  variations  was  assumed  to 
be  negligible  although  this  is  not  a  good  approximation  in  the  region  between 
the  burner  and  the  temperature  maximum.  However,  in  this  region  corrections 
tend  to  be  lower  and  limited  spatial  resolution  introduces  comparable  uncer¬ 
tainties.  For  example,  at  4  mm  in  a  0  *  3.0  flame,  the  temperature  gradient 
is  about  100  K/mm  and  the  estimated  thermocouple  correction  factor  is  200  K. 

A  50Z  change  in  the  thermal  conductivity  would  cause  an  error  in  the  corrected 
temperature  only  to  the  spatial  resolution  of  the  thermocouple.  Because  the 
calculated  ion  concentrations  derived  from  Langmuir  probe  data  are  not  greatly 
affected  by  errors  of  this  magnitude,  we  do  not  attempt  more  complicated  cor¬ 
rections  close  to  the  burner. 

Measurements  were  also  made  in  0  =■  2.25  flames  using  the  electrical  heat¬ 
ing  technique.  As  previously  shown  (Fig.  1)  a  heavily  sooted  thermocouple  (A) 
and  the  lightly  sooted  thermocouple  (B)  measurement  at  21  mm  from  the  burner 
gave  the  same  temperature  even  though  radiative  properties  (vacuum  calibra¬ 
tions)  differed  considerably.  The  value  of  €/h  extracted  for  the  lightly 
sooted  wire  (thermocouple  B)  was  used  to  correct  the  temperatures  in  two  rap¬ 
idly  recorded  unheated  thermocouple  profiles  as  described  above.  The  agree¬ 
ment  in  the  region  of  overlap  was  about  t  20  K.  Other  heating  curve  measure¬ 
ments  of  the  temperatures  nearer  the  burner,  while  somewhat  uncertain,  are  in 
accord  with  these  profiles. 

Temperature  profiles  in  leaner  flames  (e.g.,  0  =  1.5,  1.75,  2.00)  are 
complicated  by  catalytic  heating.  Use  of  either  a  heavy  Be0/Y203  coating  of 
10  to  20  pm  thickness  or  a  soot  coating  (10  pm  thickness)  eliminated  this 
effect.  Both  thermocouples  had  similar  dl/dT*  («  (€dT)°‘5)  values  measured  in 
a  vacuum.  The  assumption  of  equal  values  of  €/h  for  the  two  thermocouples 
gave  about  30-40  K  difference  between  the  two  experimental  temperature  pro¬ 
files  in  any  flame. 

Overall,  the  maximum  difference  in  temperature  measured  at  any  point  in 
any  flame  is  about  ♦  60  K  (♦  80  K  in  0  =  2.0)  while  typically  any  difference 
near  and  downstream  of  T«*»  is  closer  to  t  30  K  (±  40  K  in  0  *  2.0). 

Temperature  measurements  in  several  flames  of  interest  were  made  by  the 
above  technique  and  are  reported  in  Figs.  2  and  3.  Figure  2  gives  profiles  on 
the  AeroChem  stainless  steel  burner  at  several  equivalence  ratios,  0.  The 
flames  are  all  acetylene-oxygen  (p  ■  2.7  kPa,  unburned  gas  velocity  u  *  50 
cm/s).  In  Fig.  3  the  temperature  profile  on  the  AeroChem  stainless  steel 
burner  and  on  the  MIT  copper  burner  are  compared  for  a  0  *  3.0  flame.  In 
obtaining  these  data  three  different  thermocouple  wire  diameters  were  used, 

51,  76,  and  127  pm.  The  range  in  temperature  among  them  was  *  ±  2X.  The 
temperature  differences  between  the  two  burners,  Fig.  3,  are  certainly  due  to 
differences  in  thermal  losses,  mostly  to  the  cooled  burner  surface.  The  cal¬ 
culated  adiabatic  flame  temperature  for  these  conditions  is  2772  K.  The 
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temperature  profiles  decrease  with  different  slopes,  e.g.,  the  MIT  burner 
flame  temperature  decreases  more  rapidly  than  does  the  temperature  on  the 
AeroChem  stainless  steel  burner.  When  the  temperature  profiles  reported  here 
and  by  others’ »**  are  plotted  against  the  distance  above  the  burner  divided  by 
the  burner  diameter,  the  temperature  profiles  become  parallel.  Fig.  4. 


D.  MASS  SPECTROMETER  MASS  CALIBRATION 

Calibration  of  the  mass  spectrometer  for  mass  is  described  in  detail  in 
Ref.  22  so  will  only  be  summarized  here.  Metal  salts  were  added  to  the  flame 
to  provide  unambiguous  mass  markers  at  low  masses.  The  salts  of  potassium  (39 
and  41  amu) ,  rubidium  <85  and  87  amu),  cesium  (133  amu) ,  and  lead  (206  and  208 
amu)  were  used.  Zn  order  to  identify  the  masses  of  heavy  ions  beyond  200  amu, 
the  spectra  from  fuel  rich  deuterated  acetylene  flames  were  evaluated.  Com¬ 
parison  of  mass  spectra  from  C2H2  and  C202  f lanes  generally  showed  a  one-to- 
one  correspondence  for  the  peaks  in  the  same  general  pattern  of  relative  in¬ 
tensities;  the  deuterated  ion  peaks  were  shifted  to  higher  mass  by  a  differ¬ 
ence  that  generally  increased  for  larger  ions.  Each  of  the  ion  species  ex¬ 
hibited  a  significant  *3C  ion  peak  with  intensity,  relative  to  the  pure  »*C 
ion,  proportional  to  the  number  of  carbon  atoms,  e.g.,  1.1X  per  carbon  atom. 

Mass  assignments  were  made  to  be  consistent  with:  low  mass  hydrocarbon 
ions  identified  in  previous  work7***  alkali  metal  and  lead  ions;  and  several 
integer  mass  relationships,  AM,  between  the  ions  observed  in  the  two  isotopic 
flames.  Mass  peaks,  observed  at  every  carbon  number  from  3  to  45  (the  high 
mass  limit  of  the  instrument),  appeared  as  single  peaks  in  the  protonated 
spectra,  but  many  became  multiple  peaks  upon  deuteration. 

Our  best  estimates  from  the  mass  calibration  for  the  large  ion  identities 
and  structures  are  given  in  Table  I.  Alternate  structures  are  possible  for 
some  of  the  ions  but  given  the  C/H  ratio  and  the  structure  rules  of  organic 
chemistry  there  is  not  as  much  freedom  in  choosing  a  structure  as  one  might 
initially  suppose.  Note  that  the  odd  carbon  number  ions  can  be  represented  as 
polycyclic  molecular  ions  while  the  even  carbon  number  ions  are  represented  as 
a  protonated  polycyclic  molecule.  All  of  the  structures  are  compact  with 
little  branching.  (The  2-carbon  atom  side  chains  may  also  appear  as  condensed 
5-member  rings  in  the  polycyclic  ring  system.)  The  largest  ion  observed, 
C4sH17,  mass  557,  is  still  small  compared  to  small  soot  particles.  The  mo¬ 
lecular  ions  in  Table  I  would  most  probably  be  flat  plates  with  diameters  up 
to  about  1.0  nm  for  mass  557  compared  with  the  1.5  nm  diameter  of  the  smallest 
observed  soot  particle.3 


E.  CALIBRATION  OF  MASS  SPECTROMETER  SENSITIVITY 

Calibrations  of  the  mass  spectrometer  for  sensitivity  as  a  function  of 
mass  is  described  in  detail  in  Ref.  22  so  will  only  be  summarized  here.  One 
of  the  major  difficulties  with  mass  spectrometer  studies  of  flame  ions  is 
that  of  calibrating  the  instrument  sensitivity,  which  is  mass  dependent,  to 


obtain  ion  concentrations.  A  simplified  calibration  was  performed  based  on 
our  previous  Langmuir  probe  measurements  of  total  ion  concentrations  in  which 
we  have  great  confidence.*1 


Flame  ion  spectra,  i.e.,  the  currents  due  to  individual  ions,  were  re- 
corded  at  various  distances  above  the  burner.  Spectra  were  also  recorded  as  a 
function  of  distance  above  the  burner  while  operating  the  mass  spectrometer  as 
a  high-pass  filter  (no  applied  d.c.).  In  this  mode  of  operation  only  ions 
above  a  fixed  mass  are  sampled.  Since  there  should  be  less  discrimination 
against  high  masses  in  this  mode,  it  provided  a  better  estimate  of  the  rela¬ 
tive  concentrations  of  heavy  and  light  ions  than  using  the  individual  ion 
currents  directly.  The  integrated  currents  over  a  small  mass  range,  in  the 
high-pass  filter  mode,  divided  by  the  sum  of  resolved  ion  peak  currents  over 
the  corresponding  mass  range  was  used  as  a  measure  of  mass  discrimination. 

A  calibration  function  was  thus  derived  and  applied  to  the  normalized  individ¬ 
ual  ion  current  profiles  to  adjust  then  to  the  high-pass  filter  current 
values,  see  Fig.  5. 

The  mass  spectrometer  ion  currents  were  then  compared  with  the  Langmuir 
probe  derived  ion  concentrations  to  obtain  a  calibration.  This  was  done  by 
equating,  at  10  mm  above  the  burner  where  the  C#Hj7  ion  dominated  and  the 
curves  were  relatively  flat,  the  “sum  of  the  individual  ion  currents'*  to  the 
total  ion  concentration  measured  by  the  Langmuir  probe  in  a  similar  flame.*1 
This  gave  a  calibration  factor  by  which  all  ion  currents  (resolved  individual 
or  high-pass  filter  totals)  were  multiplied  to  give  ion  concentrations.  It 
can  be  seen  in  Fig.  6  that  both  the  current  "sums"  (of  individual  peaks)  and 
"totals"  (high-pass  filter  results)  are  in  good  agreement  throughout  the 
flame,  yet  they  deviate  strongly  from  the  Langmuir  probe  result  beyond  about 
10  mm.  The  reasons  for  this  deviation  are  discussed  in  Section  III.D. 


XXI.  RESULTS  AND  DISCUSSION 


A.  TEMPERATURE  EFFECTS 

Some  measurements  of  the  effect  of  temperature  on  soot  formation  in  pre- 
mixed  flames  were  made,  with  very  surprising  results.  This  work  has  been  pub¬ 
lished*4  and  so  it  will  be  only  cursorily  summarized  here  with  a  brief  discus¬ 
sion  of  some  of  the  implications  which  we  have  not  yet  explored.  In  these 
experiments  soot  formation  was  observed  in  atmospheric  pressure  toluene  and 
decalin  flames  of  various  fuel/0a/Na  mixtures.  As  the  0a/(0a  ♦  Nt)  ratio  was 
increased,  the  onset  of  sooting  occurred  at  higher  fuel/0t  ratios,  and  at 
higher  calculated  flame  temperatures,  consistent  with  previous  observations  of 
others.*1'**  However,  when  the  temperatures  were  measured  using  two-wave¬ 
length  emission  pyrometry,  the  flame  temperatures  at  the  soot  thresholds  were 
constant  for  each  fuel  (1750  K  for  toluene  and  1720  K  for  decalin).  These 
results  for  toluene  are  displayed  in  Fig.  7.  Soot  volume  fractions,  f«  were 
also  measured  under  a  Naval  Research  Laboratory  contract,*7  using  a  multiwave- 


length  laser  extinction  technique,  as  a  function  of  equivalence  ratio,  02/(0t 
♦  Nx>  ratio,  and  the  experimental  flame  temperature.  The  soot  volume  frac¬ 
tions  varied  strongly  with  the  0f/<0*  ♦  N#)  ratio,  but  the  variation  with 
measured  temperature  was  independent  of  the  0t/<0t  ♦  Nt)  ratio.  It,  in  fact, 
appeared  in  these  experiments  that  the  soot  yield,  for  each  fuel,  at  and  above 
soot  threshold  was  uniquely  determined  by  the  measured  flame  temperature, 
independent  of  how  this  temperature  was  obtained. 

These  results  clearly  indicate  the  hazards  involved  in  using  adiabatic 
flame  temperatures  instead  of  measured  temperatures  in  interpreting  soot  for¬ 
mation  studies.  The  results  in  Fig.  7  imply  that,  as  the  equivalence  ratio 
is  increased,  some  of  the  heat  release  reactions  do  not  go  to  completion-  - 
other  obvious  explanations  can  be  rejected.** 

Probably  related  to  the  above  is  the  equilibrium  or  nonequilibrium  nature 
of  soot  formation  in  flames.  It  is  generally  accepted  that  at  the  threshold 
for  soot  formation  the  production  of  soot  is  not  predicted  by  an  adiabatic 
equilibrium  calculation  but  that  soot  format i on  requires  a  nonequilibrium 
mechanism  for  its  formation.  Ue ,  in  fact,  use  this  as  an  argument  for  the 
ionic  mechanism  of  soot  formation  which  depends  upon  the  nonequilibrium  for¬ 
mation  of  chemi-ions  to  drive  the  process.  For  more  recent  use  of  this  con¬ 
cept  see  Refs.  28-30.  Ue  have  made  some  equilibrium  calculations  of  the  quan¬ 
tity  of  soot  formed  in  flames  at  fuel  concentrations  exceeding  the  threshold 
fuel  composition,  and  they  seem  to  predict  an  excess  of  soot  over  that  mea¬ 
sured.  This  work  needs  further  evaluation  and  substantiation  before  we  are 
willing  to  report  it  publicly}  if  substantiated  it  would  mean  completely  dif¬ 
ferent  mechanisms  might  be  involved  in  flames  near  soot  threshold  and  in  very 
rich  mixtures' 


B.  THCRWQQYNAMCS  OF  SOOT  FORMATION 

It  has  been  customary  in  the  combustion  community  to  express  the  compo¬ 
sition  of  a  fuel -oxidizer  mixture  as  an  equivalence  ratio  calculated  by  assum¬ 
ing  the  products  are  Ht0  and  C02.  The  equivalence  ratio  is  then  defined  as: 


l  .  (Fuel/Oxidizer).  «imL _ 

<Fuel/Oxidizer),»kt»i.*Ti» 

While  this  is  usually  satisfactory  it  has  been  recognized  to  be  somewhat  mis¬ 
leading  for  fuel  rich  mixtures,  such  as  encountered  in  soot  work.  Some  people 
use  the  C/0  ratio***1*  and  some  have  chosen  to  assume  the  product  is  CO  rather 
than  COf.**  In  search  of  some  guidance  on  how  to  treat  this  problem  we  have 
calculated  the  adiabatic  equilibrium  composition  at  the  soot  threshold  for  a 
set  of  fuels  spanning  the  threshold  soot  index,  TSI,  range  from  0  to  100*4'**) 
the  results  are  plotted  in  Fig.  8.  It  is  clear  that  CO  is  no  better  choice 
to  predict  sooting  than  C0t  and  that  Ht0  is  similarly  a  poor  choice. 


The  linear  variations  of  the  ratios  plotted  in  Fig.  8  are  interesting 
but  their  significance  is  not  clear}  they  are  probably  an  important  considera¬ 
tion  in  the  equilibrium/nonequilibrium  problem  poised  above! 

C.  REACTION  HEATING  OF  SOOT  PARTICLES 


In  the  ionic  mechanism  of  soot  formation  it  is  assumed  that  the  large 
ions  and  very  small  charged  soot  particles  are  formed  by  the  growth  of  smaller 
ions  via  the  addition  of  small  neutral  species,  such  as  acetylene.  Homann  has 
proposed1*  an  alternate  mechanism  to  explain  the  observed  large  ions.  He 
supports  the  hypothesis  that  as  the  incipient  soot  particles  grow,  the  ex¬ 
othermic  growth  reactions  deposit  energy  in  the  particles  more  rapidly  than 
they  can  dissipate  the  energy.  The  particle  temperature  thus  exceeds  the  gas 
temperature  so  the  particles  are  thermally  ionized.  Homann  has  not  developed 
a  quantitative  description  of  the  phenomena**}  it  is  incumbent  upon  us  to 
analyze  this  alternate  mechanism.  This  analysis  has  not  yet  been  published  so 
will  be  presented  here  in  detail. 

Ue  have  made  quantitative  estimates  of  the  amount  of  reactive  heating 
expected  for  growing  species  in  a  0  *  3.0,  2.67  kPa  C2H2/02  premixed  flame 
(50  cm  s*1  unburned  gas  velocity).  The  model  selected  is  similar  to  that  used 
by  Millikan*7  to  estimate  particle  temperatures.  Ue  first  describe  the  model 
and  quantitative  estimates  of  the  various  input  parameters,  and  then  discuss 
the  implications  on  particle  heating. 


Growth  of  species  from  molecular  size  (150  amu)  to  particle  size  (10* 
amu)  is  considered.  P, ,  designates  a  particle  which  has  undergone  n  reactive 
collisions.  The  particle  is  heated  in  the  growth  reaction  with  species  G, 


Pn-1  +  6 


kG<n> 


where  k#  <n)  is  the  reaction  rate  constant  for  the  process  creating  P. .  The 
initial  temperature  of  the  precursor,  P,  .t,  is  To,  and  the  entire  exother- 
micity,  q« ,  of  the  reaction  is  assumed  to  heat  particle  P*  of  heat  capacity, 
C.. 


Particle  cooling  occurs  either  radiatively  or  by  collisions  with  the 
flame  species,  M,  as  described  by  the  process! 


kM(n) 

P  <T.)  ♦  M  — - - >  P 

n  l  i 


-  [T*  -  -£-]  * M 


where  q«  is  the  energy  transferred  to  the  bath  gas,  and  k» (n)  represents  the 
collision  rate  constant  of  M  with  P. . 
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Several  properties  of  the  low  pressure  C2H2/02  flames  lead  to  simplifi¬ 
cation  of  the  model.  Relative  to  the  range  of  particle  sizes  considered  here 
(10  nm  diam),  the  gas  mean  free  path  is  very  large  (typically  on  the  order  of 
10  pm)  compared  to  the  particle  diameter  so  molecular  heat  transfer  processes 
dominate.  For  simplicity  we  neglect  radiative  cooling  of  the  hot  particles  so 
that  an  upper  limit  to  particle  heating  is  obtained. 


The  reactive  heating  rate  is  approximated  as  the  enthalpy  of  the  growth 
reaction,  AH*,  times  the  rate  of  reactive  collisions  of  the  growth  species  G 
with  the  particle:  AH*  k*CG3.  The  rate  constant,  k*  ,  can  be  written  as  the 
product  f*2*CG3  where  2*  is  the  collision  rate  of  G  with  the  particle  and  f* 
is  the  fraction  of  collisions  which  are  reactive.  Similarly,  the  rate  of  heat 
transfer  to  the  flame  gases  is  approximated  as  q*2NCM3  where  qN  is  the  heat 
transfer  per  collision  between  the  particle  and  M,  and  2*  is  the  collision 
rate.  qn  is  given  by: 


qH  =  a  C*  AT 


where  a  is  the  accommodation  coefficient  for  energy  transfer  from  the  particle 
to  M,  CN  is  the  heat  capacity  of  M,  and  AT  represents  temperature  difference 
between  the  particle  and  the  flame,  T*  -  T, . 


The  competition  of  these  two  processes  will  result  in  particle  heating  to 
a  steady-state  temperature  determined  by  the  condition: 


AH*  f* 2*  CGI  +  qN2N  [Ml  = 


The  particle  temperature  is  then  given  by: 


^GfGZGC6] 


V->  =  Tr 


b  CM  2  CM3 
n  n 


T* (■)  represents  the  "steady-state"  temperature  of  the  particles;  the  particle 
"temperature"  will  initially  be  a  function  of  time. 


Ue  first  consider  the  steady-state  particle  temperatures  for  various 
"particle"  sizes  from  molecular  size  (150  amu)  to  10*  amu,  corresponding  to  a 
spherical  particle  diameter,  dp,  of  6  nm  (density  s  1.5).  It  will  be  shown 
that  to  a  first  approximation,  T* (•)  is  independent  of  the  particle  size. 

Thus  the  particle  temperature  can  be  considered  independent  of  its  growth 
history  once  "steady-state"  is  reached.  The  denominator  of  the  right  hand  term 
of  Eq.  (2)  is: 


aC..2..CM3 

ri  n 


B.C.2.X. CM3 

llll 


. 
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where  the  sum  is  taken  over  the  major  flame  component s  with  mole  fractions  X,  , 
and  CM3  is  the  total  gas  concentration.  For  the  •  ■  3.0  flame,  the  major 
species  are  CO,  H*  and  C*HX  in  approximate  ratio*'*  of  6t3tl  and  CHI  *  1  x 
10* 7  molecules  cm**  at  2000  K  and  2.7  kPa.  The  collision  number  2,  from  ki¬ 
netic  theory  is: 


2. 

i 


(4) 


where  p,  is  the  reduced  mass  of  the  particle/reactant  pair.  For  large  parti¬ 
cles  this  smoothly  goes  over  into  the  rate  of  molecules  i  colliding  with  the 
surface  of  a  spherical  particle. 


2 

i 


(5) 


Even  for  a  particle  of  mass  150,  use  of  m,  rather  than  p,  results  in  less  than 
10Z  error  in  the  collision  number,  and  the  variation  in  d,  for  Hs,  CO,  and 
CXH#  (0.2-0. 4  nm)**  is  small  relative  to  d»  (•  1  nm>.  Assuming  ■  to  be  inde¬ 
pendent  of  the  gas,  we  calculate  an  effective  value  for  Eq.  (3)  at  2000  K  in 
terms  of  the  CO  diameter,  d«*  which  will  be  essentially  valid  for  all  particle 
sizes.  The  heat  capacities  used  were  estimated  to  be  (kJ  K* 1  mole**)  26,  29, 
and  75  for  Hx,  CO,  and  C*HX,  respectively.**  Likewise  the  reaction  collision 
number  2*  can  be  simplified.  Equation  (2)  is  now  (employing  m*  *  26  amu  for 
acetylene,  the  only  real  growth  candidate)40*4* 
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19  AH  f-X.(d_*  d)  CH3 
6  G  G  P  G 

«(d_  ♦  d__)*  CHI 
P  CO 


(6) 


where  AH*  is  in  kJ  mole**  and  the  diameters  are  in  angstroms.  Assuming  d*  • 
dt,  <  d» ,  or  when  d»  is  very  large,  the  size  dependence  (except  any  hidden  in 
«)  drops  out  and 
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19(AHg(k J  mole**)  f^Xg ) 
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(7) 


This  equation  holds  equally  well  for  molecules  and  particles  from  150-10*  amu 
in  a  I  ■  3.0  flame.  Now,  in  order  to  evaluate  T* ,  we  need  to  estimate  the 
heat  of  reaction  which  we  have  assumed  involves  the  major  hydrocarbon  species, 
CxHf.  We  note  that  mass  spectrometric  measurements  in  low  pressure  CfHt/0t 
flames  have  shown  that  the  next  most  abundant  hydrocarbon  species,  C«Ht,  is  an 
order  of  magnitude  lower  in  concentration.*'* 


13 


Mature  soot  particles  show  considerable  graphitic  nature  while  small, 
young  particles  have  higher  H/C  ratios  which  decrease  as  the  particle  grows. 
Ue  estimate  the  exothermicity  of  the  growth  reactions  by  two  overall  process- 


CSH«  -»  2C(s)  ♦  H, 

1/4C*H»  +  C*Ht  -»  1/4C14H10  +  1/2H* 


AH,,,  -  -227  kJ  mole*1 
AH,,,  -  -178  kJ  mole*1 


Using  available  heat  capacity  data,**»*J  at  2000  k,  the  processes  have  exo- 
thermicit ies  of  about  210  and  185  kJ  mole*1,  respectively.  Thus  Eq.  (7)  be- 


T„  S  T_  ♦  - < 

p  r  a 

for  the  more  exothermic  process. 

The  coefficients,  f,  and  a,  remain  the  major  uncertainties  as  Millikan 
found . 37 


Energy  transfer  in  reacting  molecular  systems  has  been  investigated  both 
for  thermal  and  chemically  activated  systems.  The  efficiency  of  energy  trans¬ 
fer  is  expressed  in  terms  of  an  average  energy  step  based  on  a  model  of  energy 
transfer  probability.  A  number  of  investigations  have  shown  that  the  mean 
energy  transferred  per  collision  decreases  with  increasing  temperature.**-** 
Particularly  apropos  to  the  accommodation  coefficient  a  for  present  considera¬ 
tions  are  experiments  by  Barker  and  coworkers**’*7  on  vibrational  energy 
transfer  from  highly  excited  ground  electronic  state  azulene  produced  in  in¬ 
ternal  conversion  of  the  laser  excited  molecule.  This  mimics  the  exothermic 
growth  process  in  that  the  vibrational  temperature  of  the  large  molecule  is 
much  higher  than  that  of  the  bath  gas.  Using  different  wavelengths  for  ex¬ 
citation,  the  dependence  of  the  transferred  energy  on  the  azulene  energy  con¬ 
tent  was  probed.*7  The  average  collisional  energy  transfer  to  CO  was  about 
375  cm*1  for  azulene  excitation  of  30,600  cm*1,  and  it  was  ISO  cm*1  for  17,500 
cm*1  excitation.  These  values  correspond  to  i  •  0.2  and  0.1  for  CO  in  the 
high  and  low  energy  excitations,  respectively.  The  corresponding  values 
(energy  transferred  and  a)  for  H*  (about  the  same  heat  capacity  as  CO)  were 
roughly  half  of  these.  Barker  and  Golden**  found  roughly  a  T*®* *  dependence 
for  the  average  energy  transferred  to  the  diatomic  Na  in  the  azulene  system. 

In  the  range  300-525  K  it  was  about  295  cm*1  at  room  temperature,  and  it  drop¬ 
ped  to  180  cm*1  from  400  to  525  K  in  good  agreement  with  the  observations  of 
Brown  et  al.**  for  Na  in  ethyl  acetate  from  300-800  K.  Also,  Heymann  et  al.** 
found  the  average  energy  transferred  to  CO  from  cycloheptatriene  to  exhibit  a 
roughly  T*®- *  dependence  while  other  monatomic  and  diatomics  exhibited  some¬ 
what  lower  temperature  dependencies.  Considering  the  major  post  flame  compo- 
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nent  CO,  we  correct  the  conservative  a  *  0.1  to  flame  temperatures  by  scaling 
by  the  inverse  square  root  of  the  temperature  which  gives  a  «  0.04. 

Gas-surface  accommodation  coefficients,  f • ,  are  generally  fairly  high  at 
room  temperature  (0.3-1)  and  also  have  been  observed  to  decrease  with  increas¬ 
ing  temperatures.  In  a  study  of  gas-surface  single-collision  energy  transfer, 
Rabinovitch  and  coworkers4*  found  a  vibrational  accommodation  coefficient  of 
0.2  for  cyclobutene  on  quartz  at  975  K.  We  assume  that  a  for  the  particles  is 
at  least  as  large  as  the  molecular  value  of  0.04  estimated  above. 

The  particle  growth  rate  is  obtained  from  the  electron  microscopic  stud¬ 
ies  by  Howard  and  coworkers.4  From  their  results  they  calculated  surface 
growth  rates  for  soot  particles  in  the  “standard"  acetylene-oxygen  flame  (P  * 
2.7  kPa,  u  *  50  cm/s)  below  about  2.5  pm/s.  This  corresponds  with  an  estimated 
value  of  f«  on  the  order  of  0.01  for  C*H*  as  the  growth  species.  This  com¬ 
pares  with  a  maximum  value  of  f#  ■  0.001  (for  CfHt)  from  the  atmospheric  pres¬ 
sure  flame  measurements  of  Harris  and  Weiner.40*4* 

Substituting  a  *  0.04  and  f«  *  0.01  in  Eq.  (8)  gives  a  temperature  for 
the  particle  100  K  greater  than  the  gas  temperature.  If  we  increase  the  mea¬ 
sured  gas  temperature  by  100  K  and  calculate  the  total  ion  concentration  using 
Saha’s  equation  the  results  reported  in  Fig.  9  are  obtained.  The  small  in¬ 
crease  in  equilibrium  concentration  of  charged  particles  still  gives  a  theo¬ 
retical  concentration  below  the  experimental  value.  The  most  dramatic  dif¬ 
ference  between  the  experimental  measurements  and  the  calculated  values  occurs 
at  small  distances  from  the  burner.  Here  the  measured  ion  concentration  is 
falling  while  the  calculated  value  is  rising. 

The  above  analysis  neglects  the  possible  formation  of  negative  particles 
by  electron  attachment  by  the  equilibrium 

P  ♦  e*  **  P* 

It  has  generally  been  considered  that  the  ions  are  predominantly  positively 
charged  but  both  Homann40  and  we*3  have  presented  evidence  for  negative  ions 
in  this  flame.  If  the  concentration  of  negative  ions  is  assumed  to  be  in 
equilibrium  with  the  neutral  particles  that  process  will  remove  electrons 
from  the  system  throwing  the  thermal  ionization  reaction 

P  *  P+  +  e* 

further  to  the  right,  creating  more  positive  ions.  This  has  yet  to  be  con¬ 
sidered  on  a  quantitative  basis)  we  will  do  so  before  we  publish.  The  de¬ 
crease  in  the  measured  value  with  distance  is  consistent  with  the  picture  that 
the  ions  are  formed  in  the  flame  front  by  chemi- ionization  and  then  disappear 
by  ion-electron  recombination.  At  further  distances  from  the  burner  the  mea¬ 
sured  and  calculated  charged  particle  concentrations  may  agree  because  when 
the  particles  become  larger  thermal  equilibrium  ionization  prevails.  This  is 
completely  different  from  the  ionization  we  are  concerned  about  in  the  earlier 


IS 


part  of  the  flame;  unfortunately  the  two  parts  of  the  flame  are  often  con¬ 
fused  . 

The  results  reported  in  Fig.  9  do  not  support  Homann’s  proposal.11  Me 
are  not  as  confident  as  we  would  like  to  be  however  with  the  values  of  a  and 
f«  so  have  held  up  publication  until  we  can  find  more  reliable  values. 


D.  ACETYLENE -OXYGEN  FLAMES 

1.  Total  Ion  Concentrations 

Figure  6  shows  a  great  disparity  between  the  ion  concentrations 
measured  by  the  mass  spectrometer  and  by  the  Langmuir  probe  at  distances 
greater  than  10  mm.  This  is  due  to  the  formation,  beyond  this  point  in  the 
flame,  of  large  ions  which  are  measured  by  the  Langmuir  probe  but  not  by  the 
mass  spectrometer.  The  median  ion  mass  which  we  previously  deduced*3  for 
this  flame  from  Homann  and  Stroefer’s  measurements, 50  Fig.  10,  shows  an  ex¬ 
tremely  sharp  increase,  from  55  to  A, 500  amu,  between  10  and  15  mm  distance 
above  the  burner  and  to  35,000  amu  at  25  mm  distance.  These  large  ions  are 
not  detected  by  our  mass  spectrometer.  Soot  (defined  as  particles  with  a 
diameter  greater  than  1.5  nm  which  are  detectable  by  electron  microscopy) 
first  appears  at  20  mm  above  the  burner  and  peak  in  concentration,  off  scale, 
at  35  mm.**1*  These  particles  are  located  too  far  downstream  to  be  considered 
the  precursors  of  the  relatively  small  ions  observed  here  (see  also  Ref.  1). 

2.  Ion  Spectra  and  Mechanistic  Implications 

From  the  ion  spectra  at  various  distances  above  the  flame  we  note 
that  the  ion  concentrat ions  decay  with  distance  above  the  burner,  with  the 
larger  ions  increasing  in  relative  concentration,  and  the  initial  ion,  C3H3+ 
(39  amu),  becoming  less  dominant.  The  most  common  spacing  between  adjacent 
ions  is  12  amu,  equivalent  to  the  addition  of  a  carbon  atom. 

The  occurrence  of  carbon  and  hydrogen  in  the  observed  ions  is  shown 
in  Table  II.  The  masses  of  the  ions  C*H»  +  increase  in  a  sequence  where  x  is 
an  integer  increasing  from  3  to  A5  and  where  y  is  an  odd  integer  increasing 
from  3  to  17  for  the  ions  with  the  highest  C/H  values.  The  ions  appear  to 
increase  in  mass  by  adding  simple  carbon  atoms  while  maintaining  a  fixed  num¬ 
ber  of  hydrogen  atoms  up  to  a  specific  C/H  ratio  at  which  two  hydrogens  and  a 
carbon  atom  are  added  simultaneously.  Clearly  this  sequence  does  not  repre¬ 
sent  a  series  of  elementary  chemical  reactions  but  a  sequence  of  observed 
species  which  must  be  accounted  for  by  a  set  of  elementary  reactions.  This  we 
will  do  below. 

In  Table  II,  all  of  the  ions  contain  an  odd  number  of  hydrogen 
atoms.  As  shown  in  Table  1,  the  even  carbon  number  ions  have  protonated 
structures  while  the  odd  carbon  number  ions  are  molecular  ions.  It  has  been 
argued  by  Stein*1  that  the  protonated-benzenoid-polyaromatic-hydrocarbons 
(BPAH.H+)  ions,  with  an  even  number  of  carbon  atoms,  are  in  equilibrium  with 


the  unprotonated  ion  and  a  hydrogen  atom: 


BPAH-H7  *  BPAH* T  H*  <9 

Ue  observe  the  protonated  species  which  Stein  claims  should  only  exist  at 
temperatures  below  1500  K;  the  temperatures  in  our  flame  were  greater  than 
1500  K.  Stein  attributes  our  observation  of  the  protonated  species  to  the 
rapid  rate  of  Reaction  (9)  to  the  left  in  the  cooled  flame  gases  near  the 
sampling  cone.  This  equilibrium,  if  important  (we  are  not  convinced  that  it 
is)  does  not  alter  the  arguments  presented  here  or  those  for  the  ionic  mecha¬ 
nism  of  soot  formation.  The  important  point  is  that  large  molecular  ions  are 
observed  and  their  source  must  be  accounted  for.  The  reactions  we  propose 
explain  their  production.  The  above  equilibrium  would  only  add  an  additional 
set  of  reactions  to  the  mechanism,  and  some  of  the  ion  addition  reactions 
would  involve  the  unprotonated  molecular  ion  rather  than  the  protonated  mole¬ 
cule.  Stein’s  argument,  based  on  rapid  ion  molecule  equilibria,  supports  our 
premise  that  ion-molecule  reactions  in  flames  are  very  rapid. 

Stein  did  not  treat  the  corresponding  case  for  odd-carbon  numbered 
carbon  ions,  which  we  find  to  be  unprotonated. 

The  set  of  ions  in  Table  11  can  be  used  to  define  the  types  of  reac 
tions  which  might  be  involved.  The  most  obvious  reactants  with  C3H3+  are 
acetylenes  and  our  original  reaction  schemes  included*’*?  : 


c3h3+ 

♦ 

0 

N 

X 

to 

•* 

csh3+  ♦  h2 

c3h3+ 

♦ 

-» 

C7H5+ 

c,h5+ 

♦ 

c2h2 

c,h7+ 

c,h7+ 

♦ 

c2h2 

-» 

c»  *h7+  ♦  h2 

These  types  of  reactions  account  for  the  odd  numbers  of 
removal  of  two  hydrogen  atoms;  thus  ions  of  the  general 
C„Hit.4+  where  N  is  an  odd  integer  are  observed.  As  the 
there  is  a  general  increase  in  the  C/H  ratios  up  to  2.6 


carbon  atoms  and  the 
formula  CxHh.2+, 
ion  masses  increase 
in  C*SH 


45nJ  7 


The  addition  of  polyacetylenes,  C4H2,  C6H2,  and  C2H2  (the  last  two 
reported  in  Ref.  53)  account  for  the  increased  C/H  ratio;  thus  reactions  of 
the  type 


C*H7+  + 

C4H2 

c«h2 

-> 

Ct»Hlt 

or 

C»sHn+  ♦ 

c4h2 

C:7Hj  t 

increase  the  C/H  ratio.  The  last  reaction  effectively  adds  two  carbon  atoms 
without  changing  the  number  of  hydrogen  atoms.  This  reaction  type  was  includ 
ed  in  an  earlier  proposed  ion-molecule  mechanism  of  soot  formation.**  Inter¬ 
estingly,  Warnatz  et  al.*4  in  studies  of  a  similar  flame  point  out  that  about 
35Z  of  the  C2H2  reacts  via  C4H2.  Michaud  et  al.7  deduced  equilibrium  con- 
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stants  and  heats  of  formation  of  ions  from  flame  measurements  assuming  the 
above  type  equilibria.  Bittner  and  Howard55  suggested  that  such  ion-molecule 
reactions  are  responsible  for  the  equilibrium  between  acetylene  and  diacety¬ 
lene. 


We  conclude  from  the  above  discussion  that  polyacetylenes  must  be 
considered  as  important  reactants  in  the  ion-molecule  growth  of  large  ions  in 
flames. 


We  have  yet  to  account  for  the  even  number  of  carbon  atoms  observed; 
all  of  the  reactions  above,  starting  with  C3H3+,  lead  to  an  odd  number  of 
carbon  atoms.  Reactions  of  the  type: 

kj/mol 


c3h3+  ♦  ch2=c=ch2 

- — > 

c4h3+ 

c2h4 

-  37 

- * 
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must  thus  be  important. 

Some  of  the  above  reactions  have  been  observed  (for  linear  C3H3+>  at 
near  room  temperature. 55 » 57  We  will  examine  below  whether  their  rates  in 
flames  are  sufficient  to  account  for  the  large  observed  molecular  ions. 

3.  Ion  and  Neutral  Soecies  Profiles  and  Reaction  Rate  Implications 

Individual  ion  profiles  are  presented  in  Fig.  11.  Neutral  species 
profiles  from  Bittner  and  Howard,55  Delfau  and  Vovelle,*  and  Laurendeau55  are 
plotted  in  Fig.  12  along  with  selected  ion  profiles  from  this  study.  In  gen¬ 
eral,  the  data  of  Bittner  and  Howard  and  of  Delfau  and  Vovelle  are  in  excel¬ 
lent  agreement;  where  they  disagreed,  we  report  a  mean  value,  see  Fig.  12 
capt ion. 


A  number  of  observations  are  immediately  apparent  in  Fig.  12.  Most 
striking  is  the  low  concentration  of  ions  relative  to  neutral  species.  With 
the  exception  of  C4H2,  all  of  the  neutrals  reach  their  peak  concentrations 
earlier  than  the  ion  concentration  peaks;  C4H2  peaks  at  about  the  same  posi¬ 
tion  as  ions.  The  ion  concentrations  decay  more  rapidly  than  do  the  neutrals; 
some  of  the  neutrals,  in  fact,  increase  in  concentration  with  distance.  These 
increases  in  concentration  continue  far  downstream  beyond  where  soot  particle 


formation  ceases,  thus  eliminating  reactions  of  these  polycyclic  aromatics 
with  acetylenes  as  soot  precursors2 > **  because  their  rate  of  soot  formation 
would  remain  high  beyond  where  soot  ceases  to  be  produced.  Thus,  for  example: 

Rate  -  kCC2H2HC14HB] 

would  not  change  through  the  position  in  the  flame  where  soot  is  no  longer 
formed  unless  the  slight  temperature  decrease  greatly  lowers  k.  In  contrast, 
the  coincidence  of  rapid  decay  in  ion  concentrations,  especially  larger 
ions,23  with  the  observed  increase  in  soot  concentration  is  consistent  with  an 
ion-molecule  mechanism  of  soot  formation. 

Another  striking  difference  between  neutral  species  concentrations. 
Fig.  12,  and  ions  is  the  large  decrease  in  neutral  species  concentrations  with 
increasing  mass  compared  to  that  for  ions.  For  example,  the  change  in  concen¬ 
tration  from  a  C3  neutral  to  a  Ci4  neutral  is  greater  than  three  orders  of 
magnitude.  A  similar  change  in  ion  mass  occurs  with  about  one  order  of  magni¬ 
tude  change  in  concentration;  the  ions  change  from  C3  to  C4$  ions  with  only  a 
two  order  of  magnitude  change  in  concentration.  This  might  be  interpreted  as 
indicating  a  slower  rate  of  growth  for  neutral  species  than  for  ions.  Thus, 
if  the  concentration  of  large  neutral  species  is  extrapolated  from  the  avail¬ 
able  neutral  species  concentrations,  one  would  conclude  that  large  neutral 
species  concentrations  are,  in  fact,  less  than  concentrations  of  equivalent 
sized  ions!  An  extension  of  this  type  of  argument  leads,  for  both  species,  to 
concentrations  far  less  than  the  concentration  of  soot  particles!  For  ions, 
we  know  that  the  total  ion  concentration,  when  very  large  ions  are  included, 
see  Fig.  6,  is  much  greater  than  the  concentrations  of  ions  reported  in  Figs. 
11  and  12,  and  that  the  total  number  of  ions  is  comparable  to  the  concentra¬ 
tion  of  soot  particles,  see  Ref.  23.  Equivalent  evidence  for  neutral  species 
is  not  so  definitive  (see  e.g.,  Ref.  10). 

All  of  the  ions,  with  few  exceptions,  peak  at  the  same  position  in 
the  flame.  This  is  true  for  both  even  and  odd  carbon  number  ions  independent 
of  the  suggested  structure.  This  can  occur  only  if  an  equilibrium  exists 
among  all  of  the  ions,  or  if  the  rate  of  production  of  large  ions  is  rapid 
compared  to  the  time  scale  of  the  experiment.  In  nonsooting  flames  the  ap¬ 
pearance  of  many  ions  early  in  the  flame  has  been  interpreted  as  due  to  equi¬ 
librium. so  This  again  implies  rapid  ion-molecule  reactions  at  flame  tempera¬ 
tures,  even  for  large  ions.  The  demonstration  that  such  ion-molecule  reactions 
are  rapid  in  flames  is  significant  because  there  are  few  actual  measurements 
of  ion-molecule  reaction  rate  coefficients  as  a  function  of  temperature*1 -*3 
and  to  our  knowledge  none  at  flame  temperatures. 

Ion-molecule  forward  rate  coefficients  for  nonpolar  molecules  can  be 
calculated  by  the  Langevin  equation,  which  does  not  have  a  temperature  depen¬ 
dence*1  »*3: 


where:  e  =  charge  on  ion,  a  =  polarizability  of  the  neutral  species,  and  p  = 
reduced  mass. 

Measurements  for  C3H3*  reacting  with  several  species, 56 * 67 >64  at 
near  room  temperature  are  in  agreement  with  Langevin  theory  for  the  linear 
isomer  of  C3H3+  (Table  III);  the  cyclic  structure  is  more  stable  and  reacts 
more  slowly.  These  two  structures  are  assumed  to  be  in  rapid  equilibrium  at 
flame  temperatures,  the  linear  reforming  as  it  is  removed  by  ion-molecule 
reactions.  There  is  some  evidence  in  the  measurements  of  Eyler57  that  larger 
ions  and  larger  molecules  have  lower  rate  coefficients  than  predicted  by  the 
Langevin  equation  (Table  III).  Eyler  and  associates66  have  recently  obtained 
data  which  indicate  that  the  reaction  of  C3H3+  with  C2H2  is  much  slower  at 
room  temperature  than  previously  reported.  This  does  not  represent  a  problem 
to  the  ionic  mechanism  of  soot  nucleation;  there  seems  to  be  sufficient  con¬ 
centration  of  polyacetylenes  and  other  reactants  available  to  account  for  the 
required  rates  of  ion  growth,  see  Table  III.  Arguments  have  been  made  that 
ion-molecule  reactions  are  too  slow  at  flame  temperatures  to  be  important  in 
soot  formation.  The  rapid  "equilibrium'*  among  a  large  range  of  ions  indicated 
by  the  results  in  Fig.  11  strongly  refutes  this  claim. 

The  profiles  reported  in  Fig.  11  were  used  to  estimate  the  rate 
coefficients  for  the  possible  reactions  involved  in  forming  or  destroying  the 
observed  ion.  Reactions  thus  obtained  from  the  profiles  which  give  reasonable 
rate  coefficients  with  an  assumed  neutral  reactant  partner  are  considered 
important  candidates  in  determining  the  observed  profiles.  However,  we  first 
considered  the  extent  to  which  diffusion  influences  the  profiles. 

The  rate  of  increase  of  ion  concentration  in  the  flame  front  was 
used  to  estimate  a  characteri st ic  distance,  L,  and  a  characteristic  time,  i . 
The  time  was  calculated  directly  from  the  total  gas  flow  velocity,  correcting 
for  temperature  and  molar  changes  due  to  reaction. 

The  effective  diffusion  coefficient,  D.  ,  was  then  calculated  from: 


The  characteri st ic  dimensions  of  the  experiment  are  taken  from  the  charac¬ 
teristic  dimension  of  a  cylinder66: 

w  -(i tY  y 

where  H  is  the  height  of  the  cylinder  and  R  *  radius.  H  is  taken  as  the 
thickness  of  the  profile  of  the  species  of  interest  at  half  the  maximum  con¬ 
centration;  0.4  cm  was  used  for  all  ions--certainly  within  the  accuracy  of  the 
estimate.  R  was  taken  as  the  burner  radius.  These  measured  ion  diffusion 
coefficients  from  our  experiments  are  plotted  against  ion  mass  as  points  in 


Fig.  13  and  these  are  compared  with  ion  diffusion  coefficients  calculated 
from  the  ion  mobilities,  p»  by  the  Einstein  relation*7: 


where:  k  =  Boltzmann  constant  and  the  ambipolar  diffusion  coefficient,  D.  * 

2D.  The  actual  value  is  expected  to  fall  between  these  two  curves,  depending 
upon  the  extent  of  negative  ion  formation. 

The  agreement  between  experimental  and  calculated  characteri st ic 
diffusion  times,  in  both  the  magnitude  and  shape  of  the  curve  is  surprisingly 
good.  This  tends  to  confirm  the  mass  calibration  of  the  mass  spectrometer. 

The  two  points  that  fall  farthest  outside  the  calculated  area  are  C3H3+  at 
mass  39  and  CtoH*7  <a  six-numbered  ring  with  a  four  carbon  side  chain).  The 
C3H3+  deviation  probably  relates  to  its  being  the  initial  ion.  The  reason  for 
the  deviation  of  Ct0H»+  is  not  clear. 

From  the  aaove  it  is  clear  that  diffusion  dominates  the  shape  of  the 
concentration  profiles  in  the  flame  front.  The  characteri st ic  ion-molecule 
reaction  times  must  therefore  be  small  compared  to  diffusion  times.  Recogniz¬ 
ing  this,  we  screen  candidate  reactions  involved  in  the  ion  growth  process  by 
estimating  a  minimum  value  of  the  rate  constant  for  the  candidate  reactions 
based  on  the  observed  characteri st ic  times  and  reactant  concentrations.  These 
values  are  compared  with  estimated  (Langevin  equation)  and  experimentally 
observed  rate  constants  for  some  typical  candidate  reactions  in  Table  III. 

Average  characteristic  times  decrease  with  increasing  ion  mass,  Fig. 
13:  0.55  ms  from  39-153  amu,  0.36  ms  for  165-253  amu,  0.27  ms  from  277-375 
amu,  and  0.25  ms  for  387-509  amu.  Again,  there  is  no  significant  difference 
between  odd  and  even  carbon  number  ions  in  the  same  mass  range.  Me  estimated 
required  rate  constants  (to  be  faster  than  diffusion)  based  on  a  conservative 
value  of  the  characteristic  time,  0.5  ms,  assuming  a  bimolecular  reaction  in 
which  the  neutral  reactant,  n,  is  in  excess: 


t  n 

Me  consider  two  basic  reaction  schemes  by  which  large  ions  might  be 
produced.  In  the  first  scheme  a  large  neutral  molecule  reacts  with  a  small 
ion  to  produce  a  large  ion,  e.g.: 

+  C3H3+  +  H2 

In  the  second  scheme  a  small  neutral  species  reacts  with  an  ion  to  produce  a 
larger  ion,  e.g.: 


C4H2  ♦  ^17^!+  ♦  C2H2 


These  examples,  in  fact,  represent  two  competitive  routes  used  to  explain  the 
formation  of  large  ionic  species.  Are  the  large  ions  formed  by  the  growth  of 


neutral  species  such  as  the  first  reaction,  with  subsequent  additions  of  small 
ions  to  these  large  neutrals;  or  are  the  large  ions  formed  by  the  addition  of 
small  species  such  as  acetylene  and  diacetylene  to  large  ions,  our  proposed 
mechanism  of  soot  formation  via  ion-molecule  reactions?**3*52**®  Me  are 
interested  here  in  deciding  between  these  alternate  mechanisms  for  producing 
the  large  molecular  ions  which  are  observed  in  flames. 

The  low  limits  to  the  ion  molecule  rate  constants  (calculated  from 
the  characteristic  time,  concentration  of  neutral  reactant,  and  the  half-life 
equation)  required  for  the  chemical  reactions  to  be  faster  than  the  diffusion 
process  are  presented  in  Table  III.  They  span  7  orders  of  magnitude.  For  the 
smaller  neutrals  reacting  with  C3H3+  the  agreement,  Table  III,  between  experi¬ 
ment  and  Langevin  theory  is  excellent;  for  Ct0HB  to  Ct4HB  the  experiment  gives 
a  smaller  value  than  theory.  Clearly  reactions  of  ions  with  acetylene  and 
diacetylene  are  possible  sources  of  large  ions;  only  very  modest  rate  coeffi¬ 
cients  are  required.  At  the  other  end  of  the  scale,  reactions  of  large  neu¬ 
trals  with  small  ions  are  ruled  out  by  the  unreasonably  large  rate  coeffi¬ 
cients  that  would  be  necessary.  This  occurs  because  of  the  large  change  in 
the  concentration  of  the  neutral  reactant,  n,  in  the  denominator  of  the  above 
equation  for  k  when  large  molecule  concentrations  replace  the  concentrations 
of  C2H2  and  C4H2.  One  can,  in  fact,  generalize  that  only  reactants,  in  which 
one  of  the  partners  is  in  large  supply,  can  be  effective  in  soot  producing 
reactions. 

For  reactions  of  small  ions  with  the  large  neutral  species  the  data 
in  Table  III  demonstrate  that  measured  and  calculated  reaction  rate  coeffi¬ 
cients  are  far  lower  than  those  required  to  account  for  the  experimental  ob¬ 
servations;  for  reactions  of  large  ions  with  small  neutral  species  the  mea¬ 
sured  and  calculated  reaction  rate  coefficients  are  adequate  to  account  for 
the  experimental  observations.  Thus  it  is  concluded  that  large  ions  are  not 
produced  by  addition  of  small  ions  to  large  molecules  but  rather  by  reactions 
of  large  ions  with  small  molecules.  Examination  of  the  data  for  small  neutral 
molecules  reacting  with  large  ions  reveals  that  there  are  many  small 
molecules,  including  benzene,  which  are  excellent  candidates  for  adding  to 
large  ions.  This  is  consistent  with  the  observation  of  rapidly  growing  ions. 
Thus  the  proposed  mechani sm* * 3  * 52 •**  of  ion  growth  by  addition  of  species  such 
as  C2H2  and  C4H4  to  form  larger  ions  is  confirmed.  Coupled  with  previous 
arguments, 2*3»15*23*52**2  this  is  additional  strong  evidence  that  the  mecha¬ 
nism  of  soot  nucleation  is  through  a  series  of  very  rapid  ion-molecule  reac¬ 
tions. 


E.  BEN2ENE-0XYGEN  FLAMES 

Ion  concentration  profiles  were  obtained  in  the  same  benzene/oxygen/30X 
argon  flames  in  which  Bittner  and  Howard5  measured  neutral  species  profiles. 
The  objective  was  to  increase  the  data  base  against  which  detailed  modeling 
studies  could  be  performed  and  to  compare  the  ion  profiles  in  acetylene  and 
benzene  flames.  The  same  burner  was  used  as  used  by  Bittner  and  Howard  to 
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assure  as  close  a  comparison  of  the  data  as  possible,  see  Section  II. 8.  A 
fuel  rich  but  nonsooting  flame  at  •  ■  1.8  and  a  sooting  flame  at  •  »  2.0 
were  studied.  The  total  pressure  was  2.7  kPa  and  the  total  unburned  gas  flow 
velocity  was  SO  cm/s.  These  data  have  not  been  interpreted  and  the  mass 
spectrometer  ion  profiles  have  not  been  reduced  to  absolute  concentrations. 

The  results  are  recorded  here  for  archival  purposes;  as  soon  as  time  and  funds 
are  available  they  will  be  interpreted  in  terms  of  the  mechanism  of  ion  forma¬ 
tion  in  flames  and  the  participation  of  ions  in  the  mechanism  of  soot  forma¬ 
tion. 


1 •  Total  Ion  Concentrations 

The  total  ion  concentrations  in  both  flaams  were  determined  by  Lang¬ 
muir  probe  (0.65  cm  long,  0.025  cm  diam  Pt/IOZ  Rh  wire)  measurements.  Two 
techniques  were  employed.  In  one  the  probe  current  was  measured  at  a  fixed 
negative  voltage  of  -  20  V  and  the  data  interpreted  according  to  the  theory  of 
Clements  and  Smy**.  The  other  technique  was  to  measure  complete  current/vol¬ 
tage  curves  from  -60  V  to  *18  V  and  from  these  curves  obtain  the  probe  cur¬ 
rent  at  the  plasma  potential  from  which  the  positive  ion  concentration  was 
calculated  by  Calcote’s  probe  analysis. 7»  The  currents  collected  at  -20  V  for 
these  current/voltage  curves  at  a  fixed  distance  agreed  with  the  current  pro¬ 
files  collected  at  constant  -20  V.  Complete  current/voltage  curves  were  ob¬ 
tained  with  the  hope  of  obtaining  electron  concentrations  and,  from  the  wall 
potential,  the  mass  of  the  positive  ions.  Neither  of  these  analyses  were 
successful;  possibly  more  effort  would  yield  useful  results. 

The  probe  currents  at  -20  V  and  at  the  plasma  potential,  both  ex¬ 
trapolated  from  I/V  curves,  are  plotted  in  Fig.  14  for  the  8  ■  1.8  flame  along 
with  the  positive  ion  concentration  derived  by  the  two  respective  interpreta¬ 
tions.  The  differences  in  derived  concentrations  are  greater  than  normal. 

The  higher  result  obtained  from  the  -20  V  data  is  probably  the  more  reliable, 
because  of  problems  with  this  set  of  data  in  deriving  a  plasma  potential. 
Because  we  do  not  know  the  variation  in  ion  mass  through  the  flame  we  have 
arbitrarily  assumed  a  mass  of  39  amu  in  reducing  the  data;  the  use  of  a  larger 
mass  would  increase  both  sets  of  ion  concentrations  but  would  increase  the 
lower  concentration  curve  more  than  the  higher  set  because  ion  mass  has  a 
greater  effect  on  the  interpretat ion  by  Calcote  than  the  interpretation  by 
Clements  and  Smy.  In  reducing  the  probe  data  the  flame  temperatures  measured 
by  Bittner  and  Howard1  by  radiation  corrected  thermocouples  were  used,  Fig. 

15. 


Probe  currents  measured  at  -  20  V  in  the  •  *  2.0  flame  are  reported 
in  Fig.  16.  These  will  be  subsequently  reduced  to  obtain  positive  ion 
concent r at i ons . 

Total  ion  concentrations  above  a  given  mass  were  obtained  with  the 
mass  spectrometer  by  operating  the  instrument  as  a  high  pass  mass  filter. 

This  is  done  by  turning  the  dc  voltage  on  the  quadrupole  off  so  that  all  ions 
above  a  given  mass  are  detected,  ions  below  this  mass  are  rejected.  Data  for 
the  0  *  1.8  and  0  *  2.0  flames  are  presented  in  Figs.  17  and  18,  respectively. 
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The  ion  concentration,  Measured  as  ion  current ,  by  the  Mass  spec- 
troMcter  and  by  a  Langmuir  probe,  coMpare  Figs.  17  with  14,  and  Figs.  18  with 
16,  are  very  similar  in  shape.  They  are  thus  presumably  Measuring  the  same 
ion  concentrations.  This  is  in  complete  contrast  to  the  experiments  with  an 
acet y lene/oxygen  flame,  Fig.  6,  where  the  ion  concentration  measured  with  the 
mass  spectrometer  was  much  lower  than  that  measured  with  a  Langmuir  probe. 

The  difference  was  about  300.  Why  this  great  difference  between  benzene  and 
acetylene'’  In  the  acetylene  experiment  we  argued  that  the  mass  spectrometer 
could  not  detect  large  ions  while  the  probe  could  detect  very  large  ions  and 
thus  the  difference  was  due  to  the  formation  in  the  flame  of  very  large  ions 
which  went  undetected  by  the  probe.  Application  of  this  reasoning  to  the 
benzene  flame  leads  to  the  conclusion  that  there  are  very  few  large  ions  in 
this  benzene  flame.  This  is  consistent  with  observations  by  Homann*°’7*  from 
a  completely  different  experiment. 

The  significance  of  this  difference  between  the  two  flames  for  un¬ 
derstanding  the  mechanism  of  ion  formation  m  fuel  rich  flames  and  the  role 
played  by  ions  in  soot  formation  will  have  to  await  more  detailed  analyses. 

Two  basic  questions,  which  may  be  related,  are:  why  do  the  ions  not  grow  and 
why  is  their  rate  of  recombination  different?  If  the  ions  do  not  grow,  their 
recombinat ion  coefficients  will  be  greater,  so  they  will  decay  more  rapidly. 
Or,  if  the  number  of  negative  ions  is  greater  in  the  acetylene  flame  than  in 
the  benzene  flame  the  recombination  rate  will  be  greater  in  the  benzene  flame. 
Under  such  conditions  the  ions  in  the  benzene  flame  might  disappear  by  recom¬ 
bination  before  they  have  a  chance  to  grow.  We  have  demonstrated  dramatic 
changes  in  the  recombination  coefficient  for  acetylene  flames  near  soot  incep¬ 
tion  which  substant lates  this  argument.  Homann,***  however,  reports  that  in  a 
benzene  and  acetylene  flame  producing  equal  concentrations  of  soot  the  total 
number  of  positive  and  negative  ions  in  each  flame  is  about  equal  and  the 
benzene  flame  has  about  twice  as  many  total  ions  as  does  the  acetylene  flame. 
He  also  reports  a  completely  different  distribution  of  negative  ions  with 
respect  to  mass  in  the  two  systems. 

How  the  observation  that  there  are  no  large  ions  in  sooting  ben¬ 
zene/oxygen  flames  affects  the  arguments  for  the  ionic  mechanism  of  soot  for¬ 
mation  is  not  clear.  The  rapid  disappearance  of  moderately  sized  ions  to 
form  neutral  species  which  then  grow  to  soot  particles  would  certainly  be 
reasonable.  The  main  contribution  of  the  ionic  mechanism  of  soot  formation  is 
in  the  initiation  steps  and  especially  to  account  for  the  rapid  formation  of 
polycyclic  species.  It  would  be  useful  to  have  large  neutral  species  con¬ 
centrations  in  these  two  flames  just  like  we  have  information  of  concentra¬ 
tions  for  large  ions.  More  detailed  quantitative  comparison  between  the 
acetylene  and  benzene  systems  could  pay  big  dividends  in  understanding  the 
mechanism  of  ion  formation  in  these  systems  and  the  role  of  ions  in  soot  for¬ 
mation.  The  addition  of  a  saturated  hydrocarbon  to  the  systems  studied  would 
greatly  enhance  the  value  of  this  analysis.  It  would  also  be  of  great  value 
to  have  more  quantitative  information  on  negative  ions  in  these  chemical  sys- 


The  mass  spectrometer  change  in  sensitivity  with  mass  was  determined 
by  comparing  the  change  in  ion  current  when  the  dc  was  off  with  the  sum  of  the 
peak  individual  ion  currents  over  the  same  mass  range.  These  data  were  taken 
in  the  mass  scan  mode  of  the  mass  spectrometer,  i.e.,  the  burner  to  mass  spec¬ 
trometer  sampling  cone  distance  was  fixed  and  the  rf  voltage  on  the  quadrupole 
was  varied  to  scan  a  mass  range;  this  gave  the  individual  ion  peaks  at  a  given 
distance.  The  dc  off  data  were  obtained  by  setting  the  dc  off  at  specific  rf 
voltages  corresponding  to  ions  above  a  given  mass.  A  plot  of  these  data 
against  mass  permitted  the  change  in  current  to  be  obtained  for  a  given  mass 
range.  The  results  are  plotted  in  Fig.  19.  These  corrections  have  not  yet 
been  applied  to  the  data  reported  below. 

3.  Mass  Spectra 

Ion  profiles  through  the  flame  for  the  0  *  1.8  flame  were  obtained 
by  first  sweeping  the  mass  at  a  range  of  distances  and  observing  the  ion  spec¬ 
tra  at  each  distance.  These  data  were  then  used  to  identify  the  ion  masses. 
Then  the  mass  spectrometer  was  set  on  a  specific  mass  and  the  burner  moved 
vertically  with  respect  to  the  sampling  cone  producing  individual  ion  profiles 
for  each  mass.  These  data  are  presented  in  Figs.  20-24  which  are  organized  by 
increasing  mass. 

The  masses  all  peak  at  very  nearly  the  same  distance  from  the  burner 
just  as  in  the  acetylene  flame.  The  distances  at  which  the  peaks  appear  are 
plotted  as  a  function  of  ion  mass  in  Fig.  25.  Some  of  the  masses  show  double 
peaks;  these  are  also  indicated  in  Fig.  25.  It  is  interesting  that  when 
double  peaks  appear  the  extra  peak  precedes  the  “regular"  peak.  The  mechanism 
of  formation  of  dual  peaks  remains  a  mystery.  The  maximum  ion  current  de¬ 
creases  with  increasing  mass  as  indicated  in  Fig.  26. 

Fewer  ion  profiles  were  obtained  in  the  0=2  flame  because  of  dif¬ 
ficulties  with  clogging  of  the  ion  sampling  orifice.  These  results  are  pre¬ 
sented  in  Fig.  27.  Again  dual  peaks  were  observed.  The  CjHj+  ion  does  not 
appear  until  late  in  the  flame  and  then  it  is  the  dominant  ion.  An  explana¬ 
tion  of  these  data  is  not  obvious.  The  results  are  consistent  with  previous 
results  we  reported  earlier  on  a  slightly  different  benzene/oxygen  flame.** 


IV.  PUBLICATIONS 

The  following  papers  have  been  published  based  upon  work  done  on  this 
contract  > 

1.  “Ionization  and  Soot  Formation  in  Premixed  Flames,"  D.6.  Keil,  R.J. 

Gill,  D.B.  Olson,  and  H.F.  Calcote,  Twentieth  Symposium  (International) 
on  Combustion  (The  Combustion  Institute,  Pittsburgh,  1985)  p.  1129. 


2.  "The  Effect  of  Temperature  on  Soot  Formation  in  Premixed  Flames,”  D.B. 
Olson  and  S.  Madronich,  Combust.  Flame  60,  203  (1985). 

3.  "Ion  Concentrations  in  Premixed  Acetylene-Oxygen  Flames  Near  the  Soot 
Threshold,"  D.G.  Keil,  R.J.  Gill,  D.8.  Olson,  and  H.F.  Calcote,  in  The 
Chemistry  of  Combustion  Processes.  T.M.  Sloane,  Ed.,  ACS  Symposium 
Series  249  (American  Chemical  Society,  Washington,  DC,  1984)  p.  33. 


The  following  papers  have  been  submitted  for  publication: 

1.  "Ion-Molecule  Reactions  In  Acetylene-Oxygen  Flames"  H.  F.  Calcote  and  D. 
G.  Keil,  submitted  to  Combustion  and  Flame. 

2.  "Are  Ions  Important  in  Soot  Formation?"  H.  F.  Calcote,  D.  B.  Olson  and 
0.  G.  Keil,  to  be  presented  at  the  Fuel  Division  Symposium,  Fall  ACS 
Meeting,  August  1987. 

3.  "Why  Does  Soot  Inception  Stop?”  H.  F.  Calcote,  submitted  to  Combustion 
and  Flame. 


The  following  manuscripts  are  in  preparation: 

1.  "Langmuir  and  Thermocouple  Probe  Measurements  in  Sooting  Flames,"  D.G. 
Keil,  R.J.  Gill,  D.B.  Olson,  and  H.F.  Calcote.. 

2.  "Reaction  Heating  of  Soot  Particles,"  D.G.  Keil  and  H.F.  Calcote. 

3.  Comments  on  "Effects  of  Oxygen  on  Soot  Formation  in  Methane  Diffusion 
Flames,"  by  K.  Saito,  F.A.  Williams  and  A.S.  Gordon,  Combust.  Sci.  Techn. 
47,  117  (1986),  H.F.  Calcote. 

4.  "Ionization  in  Near  Sooting  and  Sooting  Benzene -Oxygen  Flames,"  D.G.  Keil 
and  H.F.  Calcote. 


V.  PERSONNEL 


In  addition  to  the  authors  the  following  personnel  made  significant  con¬ 
tributions  to  this  program: 

Douglas  B.  Olson,  Physical  Chemist 
Robert  J.  Gill,  Physical  Chemist 
John  C.  Pickens,  Technician 
Helen  Rothschild,  Librarian  and  Technical  Editor 
(Ms.  Rothschild  also  assisted  with  data  reduction.) 

Evangeline  Stokes,  Technical  Typist 


Technical  interactions  with  other  members  of  the  scientific  community 
have  taken  several  forms,  the  foremost  being  presentation  of  our  work  at  sci¬ 
entific  meetings,  seminars,  and  workshops.  Proposals  and  manuscript  reviews, 
an  important  interaction,  require  an  increasing  effort.  Personal  contacts 
with  other  technical  people  in  this  field  have  been  maintained  by  attending 
meetings,  and  by  correspondence  and  phone  calls.  Olson  and  Calcote  visited 
Or.  Charles  Martel  at  Uright  Field  to  discuss  our  work,  especially  the  effect 
of  fuel  molecular  structure  on  jet  engine  performance. 

The  following  people  visited  AeroChem  during  this  contract  to  discuss  the 
mechanism  of  soot  formation  at  a  workshop  sponsored  by  the  Army  Research 
Office* 

T.  Brabbs  and  E.  Lezberg  (NASA/Lewis  Research  Center) 

J.  Eyler  (University  of  Florida 

U.  Flower  (Sandia  National  Laboratories) 

M.  Frenklach  and  T.  Lester  (Louisiana  State  University) 

I.  Glassman,  A.  Gomez,  and  G.  Sidebotham  (Princeton  University) 

S.  Harris  (General  Motors  Research  Laboratories) 

R.  D.  Kern  (University  of  New  Orleans) 

U.  G.  Mallard  and  R.  Santoro  (National  Bureau  of  Standards) 

D.  Mann  (U.S.  Army  Research  Office) 

J.  Tishkoff  (Air  Force  Office  of  Scientific  Research) 

The  following  presentations  were  made  based  upon  work  done  on  this  con¬ 
tract  : 

1.  "Soot  Formation  in  Combustion,"  H.F.  Calcote,  Seminar,  Stevens  Institute 
of  Technology,  Hoboken,  NJ,  9  February  1983. 

2.  "Ion  Concentrations  in  PremiMed  Acetylene  Flames,"  D.Q.  Keil.  R.J.  Gill, 
and  O.B.  Olson,  Chemistry  of  Combustion  Processes  Symposium,  18Sth 
National  Meeting,  American  Chemical  Society,  Seattle,  MA,  20-25  March 
1983. 

3.  "Ionic  Mechanisms  of  Soot  Formation  in  Flames,"  H.F.  Calcote  and  D.B. 
Olson,  AFOSR  Contractors  Meeting  on  Airbreathing  Combustion  Dynamics 
Research,  Scottsdale,  A 2,  19-23  September  1963. 

4.  "Use  of  Langmuir  Probes  in  Low  Pressure  Sooting  Flames,"  D.Q.  Keil  and 
R.J.  Gill,  Fall  Technical  Meeting,  Eastern  Section*  The  Combustion  Insti¬ 
tute,  Providence,  RI,  8-10  November  1983. 

5.  "Ionic  Structure  of  Sooting  Flames,"  H.F.  Calcote.  R.J.  Gill,  D.G.  Keil, 
and  D.B.  Olson,  Fall  Technical  Meeting,  Eastern  Section*  The  Combustion 
Institute,  Providence,  RZ,  8-10  November  1983. 
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6.  "An  Ionic  Mechanism  of  Soot  formation  in  Flames,”  H.F.  Calcote,  Invited 
Presentation,  American  Physical  Society  Meeting,  Detroit,  hi,  26-30  March 
1984. 

7.  “Ionic  Mechanism  of  Soot  Formation  -  Toward  a  Quantitative  Model, “  H.F. 
Calcote,  NBS  Workshop  on  Flame  Radiation  and  Soot,  Gaithersburg,  MO,  10 
May  1984. 

8.  “Ionic  Mechanisms  of  Soot  Formation  in  Flames,”  H.F.  Calcote,  AF0SR/0NR 
Contractors  Meeting  in  Combustion,  Carnegie-Mellon  University,  Pitts¬ 
burgh,  PA,  20-21  June  1984. 

9.  “The  Effect  of  Temperature  on  Soot  Formation  in  Premixed  Flames,”  D.B. 
Olson  and  S.  Madronich,  AeroChem  Soot  Workshop,  Princeton,  NJ,  13  July 
1984. 

10.  "Ionic  Theory  of  Soot  Formation,”  H.F. Calcote,  AeroChem  Soot  Workshop, 
Princeton,  NJ,  13  July  1984. 

11.  “Ionization  and  Soot  Formation  in  Premixed  Flames,”  D.G.  Keil,  R.J.  Gill, 
D.B.  Olson,  and  H.F.  Calcote.  Twentieth  Symposium  (International)  on 
Combustion,  University  of  Michigan,  Ann  Arbor,  MI,  12-17  August  1984. 

12.  “Soot  Formation  in  Flames”  Presented  by  H.F.  Calcote  to  six  separate 
groups  on  a  “People  To  People"  tour  through  China  in  which  Calcote  was 
the  Delegation  Leader,  9-30  October  1984. 


13.  “The  Effect  of  Temperature  on  Soot  Formation  in  Premixed  Flames”  D.B. 
Olson,  Eastern  Section:  The  Combustion  Institute,  Fall  Technical  Meeting, 
Clearwater,  FL,  December  1984. 

14.  “The  Role  of  Ions  and  Charged  Particles  in  Soot  Formation”  H.  F.  Calcote, 
Seminar  Rutgers  University,  Engineering  and  Aerospace  Dept.,  10  April 
1985 

15.  “An  Ionic  Mechanism  of  Soot  Formation  in  Flames”  H.F.  Calcote,  Invited 
Presentation  at  the  Particle  Emission  Technology  Meeting",  Naval  Post¬ 
graduate  School,  Monterey,  CA,  16-18  April  1985 

16.  “Ionic  Mechanisms  of  Soot  Formation  in  Flames,”  H.F.  Calcote,  Invited 
Speaker,  Gordon  Conference,  New  Hampshire,  July  1985. 

17.  “Are  Ions  Important  in  Soot  Formation?”  H.F.  Calcote,  Seminar,  Princeton 
University,  18  February  1986. 

18.  “Iomc  Mechanisms  of  Soot  Formation,”  H.F.  Calcote,  AF0SR  Meeting, 
Stanford  University,  California,  18-20  June  1986. 
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"Ion-Molecule  Reactions  in  Sooting  Acetylene-Oxygen  Flames,"  H.F.  Calcota 
and  0.6.  Keil,  Eastern  Section)  The  Combustion  Institute,  Fall  Technical 
Meeting,  Puerto  Rico,  15-17  December  1986. 


19. 


VII.  INVENTIONS  AND  PATENT  DISCLOSURES 


There  were  no  inventions  or  patent  disclosures  to  report  during  this 
per i od . 
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TABLE  III 

Rate  Constants  for  Reaction  of  Neutrals  with  Ions  to  Produce  Larger  Ions. 


kdO-io  cm*  s**) 
From 


Reactants _  Characteristic  Langevin,*  Experiment ally 


Neutral 

Ion* « e 

Times* 

Estimated 

Observed 

C2H2 

C3H3+ 

0.0007 

11 

10  ,12*  »' 

C2H2 

Cl4Hl 1+ 

0.0007 

8 

— 

C2H2 

C4SH17+ 

0.0007 

9 

— 

C4H2 

C3Hj+ 

0.02 

7 

10* 

C4H2 

^14^1 1+ 

0.02 

7 

— 

C4H2 

c4h3+ 

0.02 

12 

7.4* 

C3H4(allene) 

c3h4+ 

0.4 

8 

11* 

C4H4 

C4Hj+ 

0.9 

12 

11* 

CtH* 

c3h3> 

3 

16 

14 ,* 15* 

C12Ha 

CjH3* 

300 

19* 

6*  .* 

C10H« 

c  jH3+ 

400 

17 

7' 

C14H. 

c3h3-* 

1000 

19* 

6*  • 

*  Estimated  with  Langevin’s  theory  using  polarizabilities  calculated  from 
refractive  indices  and  Lorentz-Lorenz  equation. 

*  Ion  in  "estimated"  and  "observed"  reactions. 

*  Linear  C3H3+  assumed  (see  text). 

*  Independent  of  ion  (see  text). 

*  Smyth  et  al.*«. 

*  Eyler.*7 

*  Anicich  et  al.** 

*  For  1-methylnaphthalene,  CuHlQ. 

1  Eyler  and  associates**  report  a  much  lower  rate  coefficient  for  this 
react  ion. 


37 


Cv 


ii 


; 


TP-465 


0  100  200  300  400  500  600  700  800  900 

HEATING  CURRENT,  mA 

FIGURE  1  THERMOCOUPLE  HEATING  CURVES  IN  VACUUM  (OPEN  SYMBOLS) 

AND  IN  0  *  2.25  FLAMES  (CLOSED  SYMBOLS) 

Distances  from  the  burner  are  given  in  mm.  The  intersections  of  matched  solid 
curves  (see  text)  give  the  indicated  experimental  flame  temperatures.  Dashed 
curve  -  linear  extrapolation  of  thermocouple  A  vacuum  data  (see  text). 
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FIGURE  2  EXPERIMENTAL  TEMPERATURE  PROFILES  FOR  ACETYLENE- 
OXYGEN  FLAMES  FOR  DIFFERENT  EQUIVALENCE  RATIOS 

Stainless  steel  burner,  low  pressure  flames  (see  text). 
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FIGURE  3  COMPARISON  OF  TEMPERATURE  PROFILES  IN  THE  MIT 
DESIGN  COPPER  BURNER  AND  THE  AEROCHEM  STAINLESS 
STEEL  MULTITUBE  BURNER 

0  *  3.0  flames  (see  text). 
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FIGURE  4  TEMPERATURE  PROFILES  AS  FUNCTION  OF  DISTANCE 
ABOVE  BURNER  NORMALIZED  BY  BURNER  DIAMETER 

□  -  Ref.  21,  0  ■  3.5;  *  -  Ref.  9  with  diameter  from  Ref.  8  assumed; 

O  "  MIT  design  copper  burner;  ■  -  AeroChem  stainless  steel  burner. 
All  acetylene-oxygen  flames  <p  «  2.7  kPa,  u  ■  50  cm/s) .  0  -  3.0  unless 

noted  otherwise.  Diameters  indicated  in  cm. 
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FIGURE  6  COMPARISON  OF  TOTAL.  ION  CONCENTRATION  PROFILE  FROM  LANGMUIR 
PROBE  MEASUREMENTS* *  WITH  PRESENT  MASS  SPECTROMETER  MEASUREMENTS 


Adjusted  mass  spectrometer  currant  profiles  (see  text)  have  been  normalised 
by  a  constant  factor  to  make  the  sua  of  individual  ion  currents  agree 
with  Langmuir  probe  results  at  10  mm. 
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FIGURE 


COMPARISON  OF  ADIABATIC  AND  MEASURED  TEMPERATURES 
AS  A  FUNCTION  OF  SOOT  THRESHOLD,  0c ,  FOR 
ATMOSPHERIC  TOLUENE/Oa/Na  FLAMES 

-  Oa/(Oa  ♦  Na>  <  0.21J  •-  Fu*l-Air  Mixtur** 

O  —  0a/<0a  *  Na)  >  0.21. 


DISTANCE,  mm 

FIGURE  9  COMPARISON  OF  EQUILIBRIUM  CALCULATIONS 
OF  CHARGED  SOOT  (POSITIVE  ION)  CONCENTRATIONS  WITH 
MEASURED  CONCENTRATIONS 
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FIGURE  12  CONCENTRATION  PROFILES 
OF  SELECTED  NEUTRAL  AND  IONIC 
SPECIES 

OH  profiles  from  Ref.  58.  Neutral 
profiles  were  calculated  from  pub¬ 
lished  mole  fraction  profiles. 

C2H2  and  O2  profiles  are  from  Fig. 
11,  Ref.  5  (multiplied  by  2  due  to 
an  apparent  error  in  the  figure) 
and  are  indistinguishable  from 
those  of  Delfau  and  Vovelle  (DV)  in 
Ref.  9.  Other  profiles  from 
Bittner  and  Howard  (BH)  in  Ref.  5 
include  C^gHg,  C^2^a>  and  Ci4Hg. 
Profiles  of  CO,  H2,  C02,  H20,  C3H3, 
C3H4,  and  C5H5  are  from  DV.  Pro¬ 
files  for  C4H2  and  C4H4,  reported 
by  both  BH  and  DV,  are,  respective¬ 
ly  similar  in  shape  and  magnitude 
<<  factor  of  2  difference)  and  were 
averaged  for  presentation  here. 

The  corresponding  profiles  for  CSH6 
were  also  averaged  although  they 
were  significantly  different.  From 
BH,  the  CcHt  maximum  occurs  at  7  mm 
while  DV  show  the  maximum  at  4  mm. 
Also,  DV  found  lower  mole  fractions 
by  as  much  as  a  factor  of  10. 
Selected  ion  profiles  are  from 
this  work. 
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FIGURE  13  COMPARISON  OF  MEASURED  (POINTS)  WITH 
CALCULATED  (LINES)  DIFFUSION  COEFFICIENTS 
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FIGURE  14  LANGMUIR  PROBE  CURRENT  AND  ION  CONCENTRATION  PROFILES 

Experimental  currents  at  plasma  potential  (0  V)  and  at  *20  V  extracted  from 
X/V  curves  at  each  distance  in  0  *  1,8  flame.  The  calculated  ion  concentra¬ 
tions  were  based  on  theory  appropriate  to  the  collection  voltage.  See  text. 
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FIGURE  15  TEMPERATURE  PROFILE  IN  0  »  1.8  FLAME 

Experimental  temperature  profile  for  0  =  1.8  benzene/oxygen/30X  argon 
flame  corrected  for  effect  of  sampling  cone,  from  Ref. 5. 


FIGURE  16  LANGMUIR  PROBE  CURRENT  PROFILE 
Probe  currents  at  -20  V  from  I/V  curves  in  0  »  2.0  flame,  see  text 
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FIGURE  17  HIGH  PASS  FILTER  MASS  SPECTROMETER  PROFILES 

Total  mass  spectrometer  current  profiles  in  three  high  pass  filter  modes: 
passing  ALL  ions,  ions  >  110  amu,  and  ions  >  390  amu.  0  *  1.8  flame. 
Solid  curve  shifted  0.5  mm  towards  the  burner. 


FIGURE  18  HIGH  PASS  FILTER  MASS  SPECTROMETER  PROFILES 

Total  mass  spectrometer  current  profiles  in  three  high  pass  filter  modes: 
collecting  ALL  ions,  ions  >  110  amu,  and  ions  >  390  amu.  0  *  2.0.  Each  curve 
represents  an  average  of  three  separately  recorded  profiles.  The  solid  curve 
has  been  shifted  0.S  mm  towards  the  burner. 
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FIGURE  19  MASS  SPECTROMETER  SENSITIVITY  CALIBRATION  FUNCTION 

Relative  mass  spectrometer  sensitivity  as  deduced  from  ratio  of  currents 
collected  in  two  modes  of  operation  in  0  =  1.8  flame,  see  text.  Symbols 
represent  experimental  current  ratios  over  SO  amu  or  60  amu  mass  ranges 
versus  the  mass  range  midpoint.  The  curve  is  a  least-squares  cubic 

polynomial  fit  to  the  symbols. 
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FIGURE  20  INDIVIDUAL  ION  PROFILES  IN  0  -  1.8  FLAME 

Selected  mass  spectrometer  logarithmic  current  profiles  for 
individual  ions  as  indicated  by  mass  in  amu. 
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FIGURE  21  INDIVIDUAL  ION  PROFILES  IN  0  »  1.8  FLAME 

Selected  mass  spectrometer  logarithmic  current  profiles  for 
individual  ions  as  indicated  by  mass  in  amu. 
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FIGURE  22  INDIVIDUAL  ION  PROFILES  IN  0  •  1.8  FLAME 

Selected  mass  spectrometer  logarithmic  current  profiles  for 
individual  ions  as  indicated  by  mass  in  amu. 
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FIGURE  23  INDIVIDUAL  ION  PROFILES  IN  0  *  1.8  FLAME 

Selected  mass  spectrometer  logarithmic  current  profiles  for 
individual  ions  as  indicated  by  mass  in  amu. 


FIGURE  24  INDIVIDUAL  ION  PROFILES  IN  0  »  1.8  FLAME 

Selected  mass  spectrometer  logarithmic  current  profiles  for 
individual  ions  as  indicated  by  mass  in  amu. 
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FIGURE  25  LOCATION  OF  CURRENT  MAXIMA  IN  ION  PROFILES 

Distance  above  the  burner  where  individual  ion  mass  profiles  maximize 
as  a  function  of  mass  in  0  »  1.8  flame.  Species  profiles  exhibiting  two 
peaks  are  indicated  as  a  vertical  line  joining  the  two  peak  locations. 
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FIGURE  26  MAXIMUM  CURRENTS  IN  ION  PROFILES 

Maximum  mass  spectrometer  currents  in  ion  profiles  plotted  as  a 
function  of  the  ion  mass.  These  maxima  correspond  to  the  same  ions 
in  the  same  flame  as  in  Fig.  25. 


